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ABSTRACT 



A new radio direction finding system has been proposed 
which offers advantages over current systems in the areas 
of bearing accuracy and simplicity of antenna design 
requirements. This system phase codes the outputs of a 
linear array of omnidirectional antennas, and then combines 
these time-separated signals so that the sequences produced 
uniquely determine the signal arrival angle. A correlation 
process using matched filters is performed on the resulting 
sequences in order to determine the corresponding arrival 
angle . 

Resolution of the inherent ambiguities of the system, 
due to elevation angle and adjacent quadrant azimuth compo- 
nents is demonstrated using a two-array technique. A 
computer simulation of the phase coding/decoding process, 
using an Idealized signal model is used to evaluate the 
performance of the system for various input noise conditions, 
array lengths, and input signal parameters. 
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I. INTRODUCTION 



A. STATEMENT OF THE PROBLEM 

In order to determine the geographical location of a 
radio frequency-emitter. Its emissions must be detected and 
processed simultaneously by two or more radio direction 
finding (RDF) systems. The unambiguous determination of 
the angle of arrival (AOA) of RF energy received from a 
target transmitter enables each RDF system to generate a 
llne-of-bearlng along the Poyntlng vector associated with 
the transmitted energy. Ideally, the llnes-of-bearlng can 
then be plotted on a chart containing the known locations 
of the RDF sites, and the point of Intersection, or the 
area enclosed by three or more llnes-of-bearlng, will 
Indicate the approximate transmitter location. 

The accuracy with which a target transmitter can be 
located by such an operation depends upon many factors. 

The number and precision of the Individual AOA measurements, 
and the relative positions of the transmitter and RDF sites 
are probably the most Important; and each of these factors 
Is, In turn, dependent upon numerous other considerations, 
both technical and non-technlcal . The point of Intersection 
of two near-parallel llnes-of-bearlng (LOB) will vary 
significantly with relatively small standard deviation 
errors In either LOB, while the Intersection of two near- 
perpendicular llnes-of-bearlng will be Influenced only 
slightly by the same errors. It Is obvious that additional 
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AOA data, providing a wide dispersion of available lines-of- 
bearlngj can greatly enhance emitter location accuracy over 
a large geographical area. 

Increasing the number and geographical dispersion of 
radio direction finding sites has proven very effective in 
the improvement of emitter location in the HP communications 
spectrum (2 - 30 MHz). The problem of accurate AOA measure- 
ments is particularly acute in this frequency range since 
most signals of interest are received via ionospheric 
propagation. The long propagation paths involved, and the 
non-deterministic nature of ionospheric propagation, have 
necessitated the use of probabilistic techniques applied to 
multiple RDF measurements. The availability of multiple 
LOB inputs from RDF sites in a coordinated net allows a 
central processor to sort through the Inputs using an 
optimization scheme, and to choose those Inputs which yield 
a most probable emitter location with a minimum area of 
ambiguity . 

From a technical standpoint, the greater the number and 
dispersion of RDF sites, the larger the statistical data 
base from which decisions can be made. However, from 
economic and political standpoints, the number and placement 
of complex, large-array RDF systems for operation in the HF 
range are very limited. The development of a simpler, 
more compact HFDF system would allow the strategic placement 
of additional sites in those locations where more complex 
facilities do not exist, or are impractical. 
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B. STATEMENT OF PURPOSE 



The purpose of this paper is to investigate the feasibil- 
ity of a new and unique radio direction finding concept, and 
to determine the applicability of this concept to a realizable 
system for operation in the HP communications range. The 
fundamental concept was originally presented in the Masters 
thesis of R. C. Todaro [1]. The results of this original 
work indicated that the new system concept offered definite 
advantages over current RDF systems in the areas of bearing 
accuracy and simplicity of antenna design. 

The specific areas of Investigation described in this 
paper are: (1) a discussion and analysis of the ambiguities 

inherent to the system, (2) a qualitative analysis of system 
performance in the presence of background noise; and, 

(3) physical scaling of the system for operation in the HF 
frequency range. 

The remainder of this chapter contains a brief descrip- 
tion of a generalized RDF system, and a comparison of the 
RDF systems currently in use. Chapter Two contains a 
description of the new system, and a discussion of geometric 
ambiguities. Chapter Three describes the computer model of 
the system and the results obtained from the incorporation 
of background noise and scaling parameters. 

C. GENERALIZED RDF SYSTEM 

Before examining any radio direction finding system for 
its individual characteristics, it is useful to consider the 
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elements basic to any system and understand how they can 
affect its performance. A block diagram of a generalized 
RDF system is shown in Figure 1-1. 




Figure 1-1. A GENERALIZED RDF SYSTEM 
1 . Signal Source 

The signal source has only three properties which 
will affect the operation of an RDF system. These are its 
placement on the earth, its transmission frequency, and its 
type of modulation. Its location on the surface of the 
earth will determine the great circle bearing between the 
transmitter and the RDF location, but except in very simple 
cases this is insufficient information with which to predict 
the line of bearing readout from the RDF system. The trans- 
mission frequency will affect the results obtained by 
determining, in conjunction with the transmission medium, 
the mode by which the signal is sent from the source to the 
RDF site. The type of modulation will have an effect on 
those systems which compare arrival times of specific 
modulation events, or those which integrate multiple antenna 
sweeps over each line of bearing. 
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2 . Propagation Medium 



It Is possible for the propagation medium to have 
little effect on the angle of arrival of the received signal, 
but this Is not normally the case. If the signal arrives at 
the RDF site by ground v;ave propagation, then results are 
relatively simple and easily predicted. Most of the signals 
of Interest arrive via the lonoxphere, however, so Its 
effects must be taken Into account. There has been a great 
deal of work done to Investigate the character of the 
Ionosphere, and much more remains before accurate predictions 
of the transmission path of a received signal can be obtained. 
A large part of the difficulty lies In the complexity of the 
equipment necessary for measurements on the Ionosphere, and 
the large amount of time necessary to complete a set of 
measurements. Since the characteristics of the Ionosphere 
depend on the sunspot cycle as well as seasonal and diurnal 
changes, a complete set of measurements on a given path 
could not be obtained In less than eleven years. 

If accurate transmission path prediction were 
available It Is conceivable that the emitter location 
process could be performed by a single-site RDF system. 

That Is, azimuth and elevation AOA measurements could be 
used In conjunction with path prediction to compute the 
point of origin of received signals. Such a capability 
would greatly enhance the emitter location capability of a 
coordinated HFDF net. 
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3 . Signal Receptor and Encoder 



This is the first section of an RDF system over 
which some design control can be exercised, and the choice 
of this block v/ill determine to a great extent the accuracy 
and speed with which a line-of-bearing can be generated from 
the received signal. The basic purpose of this unit is to 
sense the Incoming signal, combine the proper components to 
give an output which can be uniquely interpreted as an angle 
of arrival, and present these components to a code interpreter 
for transformation to a form suitable for display. Some of 
the common types of receptors and encoders will be discussed 
in a later section describing the RDF systems currently in 
use . 

^ . Code Interpreter and Bearing Output Unit 

In most RDF systems the functions of this section 
are performed jointly by an electronic indicator sub-system 
and the human operator. Ideally, the code interpreter 
should be carefully matched to the characteristics of the 
coder so the desired output Information can be obtained in 
the shortest time and with the highest possible accuracy. 

There is, of course, no way in which the equipment associated 
with these functions can increase the fundamental accuracy 
of the system to a greater degree than is allowed by the 
components of the receptor and encoder, but it should not 
degrade their performance. An exact match occurs between 
the signal encoder section and the code interpreter of the 
proposed system through the use of a matched filter technique. 
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That is, the encoding section generates a phase coded sequence 
which is uniquely determined by the angle of arrival of the 
Incoming signal. The Inverse filter in the code interpreter 
which provides an exact match for this sequence will produce 
a unique autocorrelation (or cross correlation) output. 

This process is described in more detail in the following 
chapter. 

D. CURRENT RDF SYSTEMS 

There are four basic system types which are commonly 
used to determine the AOA of transmitted signals. The 
choice of vjhlch system to use in any particular application 
is determined primarily by the frequency range of interest, 
the degree of accuracy desired, and the type of signal 
modulation expected. These four basic types are as follows: 

1) Relative amplitude 

2) Relative phase 

■ 3) Time of Arrival (TOA) 

4) Doppler 

1. Relative Amplitude Methods 

By far the most commonly used RDF systems, these 
employ directional antennas whose patterns may be mechan- 
ically or electronically rotated. The line of bearing is 
then determined by observing the angular position of the 
antenna (or its beam pattern) when the system response is 
either at a maximum or at a minimum value. Those systems 
which employ physically small antennas (e.g., magnetic-field 
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sensitive loops) mechanically rotate the antenna while 
observing the system output. 

A loop antenna which is very much smaller than a 
wavelength has two sharp nulls in its pattern, and these 
nulls lie in a plane normal to the plane of the loop. 

Rotating an antenna of this type and noting the positions 
of the nulls will give an Indication of the bearing of the 
Incoming signal. The bidirectional ambiguity can be resolved 
by addition of the signal picked up by a separate omni- 
directional antenna located closely to the loop, but not 
so close as to affect the pattern of the loop. Addition 
of the antenna outputs will produce a cardiold pattern 
indicating the proper null to use for the bearing. 

The signal processing portion of this system (code 
Interpreter and bearing output) consists very simply of a 
standard communications receiver whose output feeds the 
deflection circuits of a circularly swept CRT display. 

Antenna position is synchronized to the CRT sweep to 
establish a common rotation rate and a fixed reference 
bearing. 

The primary disadvantage of this system is the 
relatively low overall gain of the antenna, and the fact 
that the distinct feature of the pattern (null) occurs 
where the gain is at a minimum. For these reasons the loop 
antenna is Ineffective in a low signal level/high noise 
environment. In addition, the loop antenna is subject to 
polarization errors. For ionospherlcally propagated waves. 



which have nearly random polarization, these errors may be 
as large as ±35° making the loop antenna v;orthless for long 
range DF in the high frequency band [2]. 

Instead of using a rotating element, two or more 
fixed elements may be Installed with their respective nulls 
symmetrically spaced. These individual elements may be 
loops, pairs of monopoles (thereby minimizing the polariza- 
tion error mentioned above), or any structure which will 
produce a "figure-eight" pattern with distinct nulls. The 
outputs from these elements are combined in a special motor/- 
transformei'* in which they are connected to symmetrically 
spaced stator windings. A rotating pickup coll (rotor) acts 
as the secondary winding, and its output is fed to the RDF 
receiver. The amplitude output of the i^ecelver can then 
be compared with the mechanical position of the rotor of 
the transformer, or goniometer, and the bearing of the signal 
obtained. Bidirectional ambiguities can be resolved by 
forming a cardiold pattern as in the rotating loop system. 

V/hile the figure eight patterns of the individual 
elements have relatively broad beamwidths, the combining 
action of the goniometer produces a sharper "propellor" 
pattern v/ith a beamwidth inversely proportional to the 
number of individual element pairs. This allows the use of 
the maximum gain portion of the array pattern for determin- 
ation of AOA. However, in order to achieve good figure 
eight patterns, small aperture elements (much less than a 
wavelength) must be employed. Therefore the array still 
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has a low overall gain characteristic. In addition, the 
small aperture requirement severely limits the number of 
non-interacting antenna elements v;hich may be used to com- 
prise the array. This limits the factor by which the overall 
antenna pattern beam v/ldth may be reduced, and, in turn, 
limits the angular resolving capability of the array. 

In order to improve the gain and bearing resolving 
capabilities of relative amplitude RDF systems, large 
aperture arrays have been developed. These arrays are 
usually composed of a large number of individual antennas 
located on the periphery of a symmetrical structure, ordi- 
narily a circle, and connected to a goniometer which feeds 
a single channel receiver. The code interpreter/bearing 
output section of this system is similar in design and 
operation to that associated with the simple spinning loop 
system. 

The large aperture RDF arrays that are most commonly 
in use today are of the Wullenweber design. This is a 
multiband circular array, each section of which contains 
an independent concentric array of vertical monopoles, and 
a circular reflecting screen. Two or three concentric 
sections, each with different numbers and sizes of active 
monopole elements, are designed to cover overlapping portions 
of the HF communications spectrum. The outputs from each 
of the individual antenna elements are fed to a capacitive 
coupled goniometer, which utilizes discrete rotor and stator 
pickup plates instead of Inductively linked coil windings. 
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The use of capacitive coupling, required by the large 
number of Individual antenna elements, enables the code 
interpretor/bearing output section to take full advantage 
of the fine-grain resolving capability of the antenna array. 
The goniometer forms a very narrow array beam from a set of 
40 pickup plates located on the rotor, with these 40 plates 
feeding delay lines that act to compensate for the physical 
structure of the array and transform it into an equivalent 
linear (broadside) array. 

The V/ullenweber can thus produce a very narrow 
rotatable beam providing resolving power of approximately 
3°. However, in order to achieve this resolving power, up 
to 120 separate antenna elements are required for each 
frequency range, and the physical dimensions of the array 
commonly exceed 1000 feet in diameter and 200 feet in 
height. Therefore, the array and its associated transmission 
line system is a very complex and expensive structure. 

All relative amplitude systems, regardless of array 
size, have the disadvantage that they are not Instantaneous. 
That is, the amount of time the narrow array beam spends at 
each azimuth angle is a small fraction of its total rota- 
tional period. This can produce erroneous and misleading 
results with signals of short duration such as on/off keyed 
pulse transmissions. The scan rate may be increased to 
overcome these errors, but a corresponding Increase in the 
RDP receiver bandwidth is required in order to recover the 
rapidly sampled signals. Therefore, the maximum scan rate 
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is limited by the receiver bandwidth necessary for good 
co-channel discrimination (selectivity). 

2 . Relative Phase Methods 

These systems utilize the measured signal phase 
difference between omnidirectional antennas to determine 
angle of arrival. This phase difference is given by the 
equation : 



d 



4> = 2tt-^cos 6 



(radians) 



where , 

d = antenna separation along the baseline 

X = wavelength of incident wavefront 

6 = angle of arrival, measured from the baseline. 

It is obvious that, assuming a given accuracy in 
phase measurements, the accuracy with which the AOA can be 
determined (l.e., the resolving power of the system) is 
greater for larger values of the ratio When ^ is greater 

than unity, however, there is a possible ambiguity of 2ir, 
or a multiple of 2tt, in the measured phase angle, and there 
will not be a unique solution for the AOA. 

Providing that is not too large however, the 
alternative solutions are well separated, and it is gener- 
ally possible from other considerations to select the most 
probable [3]. For example, the approximate wavelength, or 
range of wavelengths, of the incident signals can be 
determined from the tuned frequency of the RDF receivers. 
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le the separation (d) is a fixed constant. Simple logic 
cuitry can be employed to associate a range of values of 
phase shift with corresponding frequency bands of the 
receivers. 

An additional ambiguity results when the incident 
,nal wavefront has a non-zero elevation angle component 
to ionospheric propagation. The phase difference 
ation will contain an additional cosine factor of the 
vatlon AOA, and the one-equation/one-unknown situation 
1 no longer apply. The requirement for an additional 
sured parameter to uniquely determine both elevation and 
muth angles of arrival necessitates the use of a second 
n-colllnear) antenna pair. The second antenna pair is 
ally placed orthogonally to the first, resulting in a 
xlally symmetric array. Since the elevation angle 
ponent must be considered in the proposed new system, 
s topic will be discussed in more detail in Chapter II. 

The signal processing portion of the relative phase 
tern utilizes two receivers, one for each antenna pair. 

IP outputs of these receivers are feed to a sum and 
ference unit which accurately adds and subtracts the 
nals and rotates the phase of either the sum or the 
ference by Tr/2. The resultant sum/difference signals 
n drive orthogonal deflection plates of a CRT display 
t, thus achieving the effect of a goniometer rotating 
the intermediate frequency of the receivers. The display 
tern produced by this technique is an ellipse with the 
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azimuthal line of bearing along its major axis. This signal 
processing scheme produces cancellation of the effects due 
to elevation angle of arrival. 

A system of this type is difficult to build and 
maintain because of the extreme requirements on the perfor- 
mance of the receivers. Both channels must be phase and 
gain matched over the operating range in order to maintain 
the amplitude and phase relationship between the received 
signals. In addition, the frequency range that can be 
covered with a given antenna spacing is limited. The 
receiving, amplifying, and phase comparison portions of the 
system can be designed to cover a wide frequency range. 

But the antenna system cannot in practice be made to cover 
a range wider than about four to one in frequency if 
adequate resolving power without too many ambiguities is 
to be obtained [3]. 

3. Time of Arrival Methods 

The time of arrival (TOA) technique requires the 
measurement of a given modulation event at two known antenna 
locations. This modulation event may be the crest, of a 
sinusoidal modulation signal, or a specific point (leading 
or trailing edge) of a pulse signal. This technique is 
frequency independent and, thus, does not experience the 
wavelength/antenna separation ambiguities inherent to the 
relative phase technique. However, elevation angle compo- 
nents produce ambiguities Just as in the relative phase 
technique . 



20 



The TOA method is analogous to the relative phase 
method in that antenna baselines are established, and gain/- 
phase matched receivers must be used in the signal processing 
section. 

4 . Doppler Methods 

The Doppler principle of wave motion can be applied 
to a specially constructed antenna system to determine the 
angle of arrival of a signal [^]. If an omnidirectional 
antenna is rotated at a fixed rate in a circular path, the 
signal at the antenna output will be a frequency modulated 
representation of the incident signal. The peak frequency 
deviation of this signal will be: 




where 

C is the speed of propagation (meters/sec), 

R is the path radius (meters), and 
0 ) is the rotation rate (radians/sec.) 

The maximum frequency deviations will occur when the antenna 
is moving directly away (negative Af) or directly toward 
(positive Af) the emitter. Although this property is not 
affected by elevation angle components of the incident 
wavefront, the peak deviation frequency will be diminished 
by the cosine of the elevation angle and may reduce the 
resolution capability of the system. 
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is 





In practice, the antenna would not be physically 
rotated; instead a ring of omnidirectional antennas would 
be constructed and a capacitively coupled goniometer would 
sample each element. The main difference between this 
technique and the V/ullenweber system is that the Doppler 
technique samples individual low gain elements while the 
Wullenweber forms the equivalent high gain beam pattern 
of a 40-element broadside array. 

The primary disadvantage of this system is that it 
requires a circular array several wavelengths in diameter, 
thus precluding its use in the lower end of the HP frequency 
range . 
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II. A NEW RDF SYSTEM 



A. BASIC DESCRIPTION 

A new direction finding technique has been proposed 
v;hich utilizes the finite difference betv^reen arrival times 
of distinct modulation events as an incident RP plane wave- 
front propagates along a linear antenna array. Although 
the time of arrival concept is not new, the technique by 
which signals are processed, coded, and interpreted to 
extract angle of arrival information is considered to be 
unique. [1] 

Basically, the technique involves the binary phase 
coding and correlation of the time separated antenna outputs. 
The correlation of these coded signals in matched filter 
devices produces a time-compression of the signals, causing 
them to arrive at the filter output at the same point in 
time. When this occurs, a unique correlation output is 
produced. Depending on the choice of phase codes utilized, 
this unique output may be either an amplified replica of the 
original input signal (autocorrelation), or it may be a 
distinct null (crosscorrelation). 

The characteristic of this process which enables the 
determination of arrival angle is the one-to-one relation- 
ship between AOA and the particular matched filter configura- 
tion which produces this dist inquishable output. In practice, 
the outputs form a bank of matched filters, constructed to 
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match a range of arrival angles, vjould be scanned for this 
unique output. The filter producing that output would 
correspond directly to a particular arrival angle. 

The technique outlined above makes use of a special 
class of binary codes called complementary sequences. The 
time correlation of these sequences In matched filters will 
produce a unique output vjhlch Is free of the coding noise 
prevalent with all other discrete binary codes. Before 
describing In detail the application of these codes to the 
RDF system, however, a brief discussion of general phase 
coding techniques Is considered necessary for continuity. 
Therefore, these techniques, and the specific advantages of 
complementary sequences, are presented In the following 
section . 

B. PHASE CODING TECHNIQUES 

Phase coding techniques were first used (In communica- 
tions systems) to Improve the range capability of pulse 
radars. V/hlle not relevant to the topic of this paper, 
this utilization of phase coding provides a graphic descrip- 
tion of general coding techniques and an example of their 
application to an actual system. 

Two parameters of a transmitted radar signal may be 
varied In order to Increase the total RF energy Illuminating 
a radar target. These are the peak pulse power, and the 
pulse duration. The peak power of the transmitted pulses 
may be Increased by utilizing higher power amplifiers In 



the transmitter output stage. However, this has proven 
inefficient from a power standpoint, since many radar 
transmitter tubes are peak power limited vjhile operating at 
only a small fraction of their average power capabilities. 

The average power output of these tubes may be increased by 
widening the transmitted pulses. However, this too has a 
drawback because the range resolution capability of a pulse 
radar is inversely proportional to pulse width. 

As a result of this dilemma, intrapulse modulation 
techniques were developed so that pulse widths could be 
increased without sacrificing range resolution. The most 
common example of intrapulse modulation is the "chirping" 
or pulse compression system in which the RF carrier frequency 
of each pulse is swept between fixed limits during its trans- 
mission. The received pulses are then passed through a 
matched filter whose response characteristic (delay versus 
frequency) is the exact inverse of the transmitted pulse 
characteristic (frequency versus time). This allows all the 
frequency components of the wide reflected pulse to arrive 
at the output simultaneously, resulting in a much narrower, 
amplified pulse signal to the detector stage. Phase coding 
of the transmitted pulse is another technique of intra- 
pulse modulation. 

The phase coding technique utilizes a discrete process 
to encode pulses so that each return pulse, when input to a 
matched filter, causes discrete, coded segments of the pulse 
to arrive at the output simultaneously. Each radar pulse is 



25 



divided into "n" equal time segments, and a coded sequence 
of "phase states" is assigned to these segments. The states 
determine the phase relationship between the RP carrier 
during that segment, and a fixed phase carrier reference. 

If the only allowable phase states in each segment are 0° or 
180 *^, then the coding is termed binary phase coding. 

Although multiphase coding schemes have been devised and 
used in certain applications, only binary phase coding will 
be discussed here. 

A representation of a binary phase coded pulse is shown 
in Figure 2-1 (a). The RP carrier is divided into seven 
equal time segments, v;lth each segment assigned a plus ( + ) 
or minus (-) phase state. The phase state values are here 
defined as no carrier phase reversal (+) and a l80° carrier 
phase reversal (-). A unit amplitude is assumed for the 
pulse . 

The filter which provides an exact match for this coded 
pulse is shown in Figure 2-1 (b), v;here the discrete delay 
units are equal to the pulse segment durations and the (+) 
and (-) phase reversals occur in reverse order from the 
transmitted sequence phase states. In a physically realizable 
filter, the discrete time delays would be produced by delay 
lines or electronic delay circuits, and the phase control 
would be provided by buffer amplifiers with inverting or 
non-inverting outputs. For this example, the buffer/phase 
control devices are assumed to have unity voltage gains. 
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Figure 2-2 is a time diagram showing the output waveforms 
at the intermediate points "a" through "g", and the final 
summation output waveform. 

There are two characteristics to be noted about this 
output. First, it contains a narrov/ main pulse that is 
delayed n-1 time units from the leading edge of the received 
pulse. Since this narrow, main pulse would be utilized by 
the radar receiver's detector and timing circuits to deter- 
mine target range, this delay must be compensated for. 
Second, there are undeslred secondary pulses in the output 
waveform which might appear to the radar receiver as excess 
noise or false targets. 

There are other codes, called pseudo-random noise codes, 
which will produce smaller secondary pulses than those 
produced by the (randomly selected) code utilized in this 
example. However, all single-sequence coding schemes of 
this type produce a certain amount of "coding noise", as 
these secondary pulses are called, and therefore exhibit a 
finite processing signal-to-noise ratio. 

There is a method, however, which will completely 
eliminate the coding noise, although it requires the simul- 
taneous transmission of two distinct coded sequences. These 
codes, called complementary sequences, v;ere first utilized 
by M. J. E. Golay in the field of multislit spectrometry. 
Since then, they have been utilized in various radar and 
sonar system designs. Golay mentioned in 196I that these 
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sequences might be applied to data transmission [5], and 
there has been some work toward that end [6]. 

A great deal of work has gone into the investigation of 
properties of complementary sequences, and the methods by 
which they may be generated. However, only the basic 
characteristics of these sequences will be presented here. 
For more detailed studies, the reader is referred to 
References 5 and 7 listed at the end of this paper. 

A complementary sequence pair is defined to be a pair 
of equally long binary sequences which have the property 
that the number of pairs of like digits (or elements) with 
any given separation in either sequence is equal to the 
number of pairs of unlike elements with the same separation 
in the other sequence. 

For example, considering the following sequences 

A: _ 

B: _-_+++_+ 

the number of like elements in A with separation distances 
of one element is three: 

A: _ 

lluuluu 3 like elements 

The number of unlike elements in B with separation distances 
of one element is also three: 

B: ___+++_+ 

llulluu 3 unlike elements 
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Similarly, the number of like elements in A and the number 
of unlike elements in B for all separation distances are 
listed belovt: 

Separation Distance # Likes in A § Unlikes in B 

2 3 3 

3 ^ ^ 

2 2 

5 2 2 

6 11 

7 11 

Golay has shown [5] that the sequence lengths must 
always be an even number, and must be the sum of tvjo squares. 
For example, the possible sequence lengths up to 50 are: 

2, 3, 8, 10, 16, 18, 20, 26, 3^, 36, 40, and 50. He has 
also shown [5] that 2^ complementary sequence pairs (some 
identical) can be generated by performing the following basic 
operations on a given sequence pair: 

1. Interchanging the sequences 

2. Reversing the order of the first sequence 

3. Reversing the order of the second sequence 
Complementing the elements of the first sequence 

5. Complementing the elements of the second sequence 

6. Complementing the even-order elements of both sequences 

An example of a pair of complementary sequence coded 
pulses being processed by a matched filter pair is shown 
in Figures 2-3 and 2-H . The input pulses, utilizing comple- 
mentary sequence codes of length eight, are shown in Figure 
2-3(a). The pulses and their matched filters, shown in 
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Figure 2-3 (b), are labeled ’A' and 'B' respectively. The 
outputs from these matched filters, and the final combined 
output, are shown in Figure 2-H. This final output waveform 
shows the signal processing advantage of this coding scheme. 
The total energy content of both received pulses has been 
compressed and combined into a single narrov/ pulse with an 
amplitude gain factor of twice the code length. In addition 
the complete cancellation of all undesired coding noise 
indicates that an infinite processing slgnal-to-noise ratio 
can be obtained with this technique. 

In radar system applications, this property allows trans 
mltters to be operated at high average power levels while 
utilizing pulses with high frequency content. Thus, both 
transmitted energy and range resolution goals may be met. 

The price that must be paid for this advantage is the 
increased complexity of the signal generation, transmission, 
and detection components. 

The signal coding/decoding process illustrated above may 
also be described in mathematical terms. If we let a^ and 
be the elements of any pair of n-element binary sequences 

f ! 

(taking on values +1 or -1), and let a^ and b^ be the phase 
reversal portions of their corresponding filters, then the 
output sequences from the filters can be represented by: 



n-j , 

c . = E a . a . , . 
J 1=1 i i+J 



n-J , 
" ifl 
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where the subscript "j" takes on integer values from -n+1 
to n-1 and, by definition: 



c 



-J 





= d 



J 



The summation equations represent a discrete convolution 
of the phase codes and their respective filter codes, 
resulting in correlation of the two sequences. In the case 
where the filters have "matched filter" characteristics 
(i.e., a^ = a^^) then the sequences are said to be auto- 
correlated. These statements apply to any pair of equal- 
length binary sequences. Then if the input sequences 
utilized are a complementary coded pair, and matched filters 
are used, the following conditions apply: 



s 



J 



c . 

J 




for j 5 ^ 0 



s = c + d = 2n 
o o o 



for j = 0 



This is a mathematical description of the final output 
sequence obtained in the process illustrated above. The 
final summation terms, s., are zero for all values except 

i) 

j = 0, which corresponds to the main pulse output during the 
eighth time segment. 

A second useful property of complementary sequences is 
that, for each sequence pair, there is an "orthogonal" pair 
which produces a unique crosscorrelation output. That is. 
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1 ^ 










the output sequences c. and d. obtained by correlating a 
complementary sequence pair with its corresponding orthogonal 
pair will be described by the relation: 



Sj = Cj + dj = 0 for all values of j . 

This process would not be useful for pulse radar or data 
transmission systems, since the processed signal has zero 
amplitude for all values of time. Hovjever, the cross- 
correlation output is applicable to the RDP system, since it 
provides a unique correspondence between the angle of arrival 
and the signal processor output. 

The generation of an orthogonal sequence from a given 
complementary sequence utilizes four of the six basic 
operations described earlier in this section. The generation 
algorithm, shown below for an eight-element complementary 
sequence pair, is given by French [6] in his Master’s thesis. 

(1) Given a pair of complementary sequences: 

A= + _ 

B= — — — + + + — + 

(2) Perform operation §1 (interchange sequences): 

A' =---+++-+ 

B' = + + - 

(3) Perform operations §2 and §3 (reverse the order of 
the sequences): 

A” = + — + + + — — — 

B”=_+ +___ 
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(^1) Perform operation (complement elements of A”): 

A'” = - + ___ + + + 

B”’ 

The resultant sequence pair A'*', B' ' ' is orthogonal to the 
original sequence pair A,B. 

An example of the crosscorrelation process, utilizing 
orthogonal complementary sequence pairs, is shown in Figures 
2-5 and 2-6. The phase coded pulses in Figure 2-5 (a) 
utilize the eight element complementary sequences A and B. 
The correlation filters, shown in Figure 2-5 (b), utilize 
the orthogonal sequences »**a**' and **’B***. The outputs 
from these filters, as shown in Figure 2-6, add to zero for 
all values of time. 

With this introduction to complementary sequence coding 
and an application to signal processing, the basis is 
established for its application to the proposed RDF system. 

C. APPLICATION OF COMPLEMENTARY SEQUENCE CODING TO AN 
IDEALIZED RDF SYSTEM 

As previously mentioned, the phase coding method of 
determining the angle of arrival depends upon the time 
difference between antenna responses in a linear array. 

This time difference is a unique function of the inter- 
antenna spacing, and the arrival angle of the wavefront. 

This is illustrated in Figure 2-7(a), in which a single 
rectangular RF pulse is assumed Incident to a four element 
array . 
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The rectangular RP voltage pulses induced at the antenna 
outputs are separated in the time domain by a shift (At) 
given by the equation: 

At = — cos 6 
c 



where , 

d is the antenna spacing, 

c is the propagation velocity of the wavefront, 
and 6 is the arrival angle of the Poynting vector 

associated with the wavefront. 

Figure 2-7 (b) illustrates the relative time spacing of the 
induced RP voltages, with the output of antenna number one 
chosen as a reference. For a signal incident upon the 
antenna system at an arrival angle of zero (along the array 
baseline) the time difference will be equal to the maximum 
propagation time betv/een antennas. The time difference 
will reduce to zero as the arrival angle Increases to the 
point where it is normal to the array baseline. 

The type of antenna used in the array is not a critical 
factor in the RDF operation of the system; any omnidirectional 
antenna can be used. Hov:ever, the antennas should be identi- 
cal in their electrical characteristics (pattern, gain, and 
Impedance), and they should be sufficiently separated to 
minimize mutual coupling. 

The first step in the phase coding technique is the 
transformation of these time shifted signals into a pair 
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of complementary sequences that can be interpreted by 
matched filter decoders. The basic hardware necessary for 
this process is shown in Figure 2-8. Each antenna feeds 
a dual output multicoupler device which can provide both a 
normal phase output and an inverted phase output (relative 
to the input signal). The outputs from these devices are 
configured so that their phase states form a complementary 
sequence pair. For this reason, the array length must 
correspond to one of the possible complementary sequence 
code lengths (2, 8, l6, l8, 20, etc.). The four element 

array was chosen here for simplicity. 

After phase coding, the output signals corresponding 
to each of the complementary sequences are combined in two 
signal channels and applied to a code interpreter. The 
code interpreter consists of a pair of correlation filters 
which, when perfectly matched to these input signals, 
produce unique autocorrelation (or crosscorrelation) outputs 

To illustrate the phase coding and decoding processes 
in more detail, the four element array is used with an 
incident wavefront at an arrival angle of zero degrees. 

The complementary sequences are identical to those in 
Figure 2-8, and an idealized rectangular RF pulse is assumed 
with a pulse width equal to twice the inter-antenna propa- 
gation delay (2—). To simplify the present example, it is 
also assumed that the phase reversal and summing operations 
are performed at the carrier frequency of the incoming 
pulse signal, and that no other signals are present. The 
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practical requirements for signal preselection, and Inter- 
slgnal discrimination circuits are considered In a later 
section of this chapter. 

Based on the above conditions, the signals at the 
antenna outputs v\^ould be as shovm In Figure 2-9. The output 
of antenna number one Is again chosen as a reference, and 
the A & B outputs of each antenna correspond to the channels 
Into which the phase-coded outputs are combined. For 
convenience, the peak carrier amplitudes of the wave forms 
are normalized to unity, and upward or downward pointing 
arrows are used to Indicate their respective phase states. 

An upward pointing arrow Indicates a non-lnverted output, 
and a downward pointing arrow Indicates an Inverted output. 
When these Individual antenna outputs are grouped and then 
added In an Ideal linear (and lossless) summer, the resultant 
output waveforms are assumed to be as shown In Figure 2-10. 

It should be noted that the amplitudes and phases of 
these Idealized waveforms are based on the assumption that 
the antenna outputs are combined with perfect carrier phase 
cancellation (l80° phase difference) or perfect phase rein- 
forcement (0° phase difference). To assume that perfect 
phase matching occurs for all possible values of time differ- 
ence (and, therefore, for all possible values of relative 
phase difference) Is a questionable assumption. That Is, 
perfect l80° or 0° matching between phase coded antenna 
outputs will occur only for those arrival angles which 
produce phase differences that are multiples of 2ir. 
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However, the waveforms in Figure 2-10 do provide a 
simplified Illustration of the basic phase coding/decoding 
technique applied to the RDF system. In addition, it will 
be shown in the follov/ing section that the relative carrier 
phase differences between individual antenna outputs do 
not affect the validity of the overall process so long as 
the signals are combined in a linear device and the principle 
of superposition can be applied. 

Returning then to this simplified example, the coded 

and summed signals occupying channel A and channel B are 

each applied to correlation filters which are configured 

as shown in Figure 2-11. The filters consist of discrete 

delay devices (or tapped delay lines) and dual output 

buffer devices with inverted and non-inverted outputs. The 

dual outputs of the buffer units are configured so that both 

autocorrelation and crosscorrelation of the coded input 

signals can be performed. For perfect f ilter-to-signal 

matching, the filter time delays must equal the inter- 

antenna propagation time (— cos 0). For this example, 

c 

(e = 0°) perfect matching requires that the time delays 
d 

equal — . The autocorrelation and crosscorrelation filter 
c 

outputs are shown in Figure 2-12 for the perfect matching 
case. As in Figure 2-9, upward and downward pointing arrows 
indicate the signal phase states, and the amplitudes corres- 
pond to peak values. When the corresponding (autocorrela- 
tion and crosscorrelation) outputs of channels A and B are 
combined in a linear device, the final output waveforms will 
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be as shown In Figure 2-13. The autocorrelation output, 
for perfect filter-to-signal matching, is an exact replica 
of the original RF pulse, with a processing gain equal to 
twice the code length. The crosscorrelation output is 
everywhere zero because of the complete phase cancellation 
of all signal components. 

A computer model of this process has been used [1] to 
determine the effect of correlation filter mismatch upon 
the autocorrelation and crosscorrelation output waveforms. 

An important result of that study was the determination 
that the crosscorrelation output is uniquely zero for the 
perfectly matched case, and that a very slight amount of 
mismatch will cause a rapid Increase in the crosscorrelation 
output. It was concluded, then, that the distinct null 
effect produced by the crosscorrelation process could be 
utilized to accurately determine the angle of arrival of 
incident signals. 

The computer model utilized in that study employed the 
idealized signal model that was assumed above, and it also 
assumed a noise-free signal environment. In order to more 
fully evaluate the performance of the complementary sequence 
phase coding process for this paper, the computer model was 
modified to Investigate the effects of background noise and 
various signal parameters. In addition, the model was 
utilized to investigate the effects of array length upon the 
overall process. A description of the model, and the results 
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and conclusions reached, are presented in Chapter III of 
this paper. 

D. SINUSOIDAL ANALYSIS 

The following sinusoidal analysis of the coding/decoding 
process will show that the basic technique schematically 
described in the previous section is valid for all arrival 
angles (between 0° and 90°), and that a 0° or l80° phase 
relationship between the antenna outputs is not a require- 
ment for proper operation of the system. 

First of all, it is assumed that the physical constraints 
on the system (array size, phase coding sequences, matched 
filter configuration, and signal parameters) are identical 
to those specified in the previous example. However, for 
the purpose of this analysis, the arrival angle may take on 
any arbitrary value between 0° and 90°. 

Therefore, the time differences between antenna responses 
will be 

At = — cos 6 
c 

and the relative phase differences will be 

A(j) = 2ir Y cos 0 = w ^ cos 6 , 

A C 

where w is the carrier frequency of the Incident signal. 

Now, if antenna number one is chosen as a time and phase 
reference, then the signal voltages induced in the four 
antennas may be represented by the equations : 
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d 

V-J^(t) = [u(t )-u(t-2^) ]cos wt = u^(t)cos wt 

V 2 (t) = [u(t-At )-u(t-2^-At ) ]cos (o3t-A(})) = U 2 ( t ) cos ((jjt-A(}) ) 

v^(t) = [u(t-2At )-u (t-2— - 2At ) ]cos ( wt-2A<})) = u^(t )cos(wt-2A4)) 
0 C D 

V 2 |(t) = [u(t-3At)-u(t-2^- 3At )]cos(wt-3A<l)) = ( t ) cos (wt-3A4> ) 

These equations represent the signal waveforms shown in 
Figure 2-7. The gating functions u^(t) through U|^(t) define 
the leading and trailing edges of the signal pulses, which 
are of duration 2 — . These signals are then phase coded 
using complementary sequences and combined in a pair of 



linear 


devices 


so 


that 


the 


resultant outputs are represented 


by : 














Channel 


A: 






+ Vj + Vj - 




Channel 


B: 




1 — 1 
> 
+ 

II 


+ Vj - V3 + V, 



The composite signals and Vg are then input to auto- 
correlation and crosscorrelation filters with discrete time 
delay segments equal to — cos 0. Since the time delays are 
equivalent to negative phase shifts (-A<f)) equal to cos 6), 

then the four correlation outputs can be represented by the 
series: 
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Channel A (autocorrelated) : 



Channel B (autocorrelated): 

~ "^^B “ '*' + v^ /-3A(}) 



Channel A (crosscorrelated) : 

\ ^ ""^A " ■*■ “ v^A3M 

Channel B (crosscorrelated): 

Vg = +Vg + Vg /-A({) + Vg /-2A({) - Vg /-3A(j) 



where v^ /-nAtj) represents the composite voltage (v^ or Vg) 
phase shifted by an angle -nA<j). For example, expanding 
v^ /-nA(j) produces the follovjing equation: 

v^ /-nA(f) = Vj^ /-nA(}) + v^ /-nA(j) + v^ /-nA(j) - V|| /-nA<{) 

A general expression for the series terms is 



Vj, /-nA(j) 



= { 



Uk(t) 



cos[wt - (k-l)A(j)] 



I /-nAtj) 



(k=l,2,3,^) 

(n=0,l,2,3) 



where u^(t) = u[t - (k-l)At] - u[t - 2 (k-l)At] . 



Since the phase delay /-nA({) is equivalent to the time delay 
nAt , the gating function u^(t) is delayed by nAt , and becomes 



Uj^(t+nAt) = u[t - (k+n-l)At] - u[t - 2 ^ - (k+n-1 ) At ] = 



ill 



Therefore, the general series term may be expressed as: 
Vj^ /-n A(}) = cos [w t - (k+n-1) A(}) ] 

Now, applying these terms to the autocorrelation and cross- 
correlation series: 



(k=l) (k=2) (k=3) (k=i^) 

= -U^COS (ojt ) - U2COS (d)t-A(j) ) - U^COS (o)t-2A(i) ) + u^^cos (uit-3^<P ) 

+U 2 COS (wt-A(J)) + u^cos (o)t-2A(J)) + U|^cos (o)t-3A<}> ) - u^cos (o)t-4A(t) ) 
+U2Cos(wt-2A(|)) +U2jCos(o)t-3A(})) +u^cos(o)t-4A(})) - UgCos(oL)t-5A(|)) 
-U2|COs(wt-3A<l)) +u^cos(o)t-^A(l)) +UgCOs(a)t-5A4)) -u,^cos(a)t-6A<t)) 



(n=0) 

(n=l) 

(n=2) 

(n=3) 



= -u^cos(iot) + u^cos (wt-2A4> ) + 4u|jCos(wt-3A((>) +u^cos(o)t-4A4)) -Uycos(o)t-6A(})) 



(k=l) (k=2) (k=3) (k=4) 

Vg = +u^cos (mt ) +U 2 cos(o)t-A(t)) -U 2 COs(wt- 2 A(()) +U2jCos(wt-3A(t)) (n=0) 

-U 2 COS (wt-A(})) - U2Cos(wt-2A(})) +U|^cos(o3t-3A(j)) - u^cos(o)t-4A(})) (n=l) 

+U2COs(wt-2A(J)) + U|^cos(wt-3A(j)) - u^cos((ot-4A4)) +UgCOs(mt-5Ai})) (n=2) 

+U|jCos(a)t-3A(})) + u^cos(wt— 4 a<})) -Ugcos(o)t-5A(})) + UyCos(mt-6A4)) (n=3) 

= u^cos(a)t) -u^cos(wt-2A(J)) + Au|^cos(mt-3A4>) -u^cos(wt-^A(})) +u,^cos(wt-6A(()) 
(k=l) (k=2) (k=3) (k=^) 

f 

= -u^cos(wt) 

-U2Cos(wt-A4>) -u^cos(a)t-2A(})) -U|jCos(oJt-3A<}>) +u^cos(wt-4A(})) (n=l) 

+U2Cos(wt-2A(})) +U£|Cos(wt-3A<l>) +u^cos(wt-^A(})) - ugCos(wt-5A4>) (n=2) 

-U|jCos(wt-3A(()) -u^cos(wt-^A(l)) - ugCos(wt-5A4>) +u^cos(ti)t-6A(l)) (n=3) 

= -u^cos(wt)-2u2Cos(cot-A(l))-U2Cos(wt-2A(f))+u^cos(mt-i}A(}))-2ugCOs(wt-5A(J)) 
+UyCOs(wt-6A4)) . 



U 2 COs(wt-A(l)) - U 2 Cos(wt- 2 A(})) + u2jCos(mt-3A4>) (n=0) 



(k=l) 



(k=2) 



(k=3) 



(k=iO 



t 

Vg = +u^cos(wt) +U 2 cos(cot-A(|)) -U2cos(cot-2A({)) +U|^cos(wt-3A({)) (n=0) 

+U 2 Cos(wt-A({)) + U 2 cos(cot- 2 A 4 >) - U||COs((ot-3A({)) + u^cos(cot-ilA4)) (n=l) 
+U2Cos(wt-2A({)) +Uijcos(cot--3A({)) -u^cos(wt-^A({)) +Ugcos(tot-5A({)) (n=2) 

-U|^cos(wt-3A(})) - u^cos(cot-^A({)) +UgCos(wt-5A(J)) - UyCos(wt-2A(t)) (n=3) 

= u^cos(wt) + 2u2COs(wt-A4)) + U2Cos(cot-2Acf)) - u^cos(ut-^Acf)) 

+ 2ugCos(wt-5A(|)) - UyCos(wt-6A(j)) . 

Finally, when each channel's autocorrelation and cross- 
correlation terms are combined, the resulting outputs will 
be : 

Vo = + Vg = 8ui^cos(wt-3A4>) 

= 8|u[t-3At] - u[t - 2 ^-3At]|cos(wt-3A(f)) 

v„ = v; + v; = 0 

These equations represent the final output waveforms shown 
schematically in Figure 2-13. The autocorrelation output is 
an exact replica of the Incoming pulse signal, and is 
amplified by a factor of twice the code length. The gating 
function Uj^(t) corresponds to the result in Figure 2-13 
in that the leading edge of the output pulse is delayed 
3At time units from the initial intercept of the incident 
pulse. Also, as in the graphical example, the final 
crosscorrelation output is everywhere zero. Therefore, it 
is concluded that the phase coding process is valid for all 
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arrival angles between 0° and 90°, that it provides a signal 
processing gain (in the autocorrelation mode) of twice the 
chosen code length, and that it provides a crosscorrelation 
output that is everywhere zero. 

E. ENGINEERING CONSIDERATIONS 

In order to apply the basic phase coding technique 
described above to a physically realizable RDF system, 
certain practical problem areas must be considered. During 
the investigation and evaluation of this technique, it 
became obvious that several factors have a direct effect 
upon the operation of the proposed system. Specifically, 
those factors that were considered most important are: 

(1) possible sources of error resulting from angular 
ambiguities, (2) the requirement for signal selection and 
discrimination, and (3) the effect of different modulation 
types on the correlation processes. The constraints and 
physical requirements which these factors place upon the 
design of such a system will, in large part, determine the 
feasibility of utilizing this type of RDF system over some 
other type. 

The first two factors (angular ambiguity and signal 
preselection) are discussed and evaluated in this section. 
The actual design requirements that must be imposed to 
satisfy these two areas are straightf orvrard , and are 
considered to be technically feasible. However, the overall 
practicality or economic feasibility of these requirements 



will depend, in part, upon the particular application of 
the system, and are not addressed in this paper. Modulation 
effects are best described and evaluated by the computer 
model used to simulate the phase coding and decoding portion 
of the system. This discussion is presented in the 
following chapter. 

1 . Angular Ambiguities 

Until now, the determination of arrival angle has 

been assumed to be a two-dimensional problem, in which one 

single-variable equation is evaluated to yield a single 

valued unambiguous solution. However, in most environmental 

situations, the direction of arrival of an Incident wave- 

front will generally be composed of both azimuthal and 

elevation angle of arrival components. The additional 

consideration of elevation angle Introduces a third spatial 

dimension to the previous two dimensional problem. That is, 

the elevation angle introduces a second variable to the 

time delay equation, so that At = — cos 0 cos a, where a 

c 

is the elevation angle of arrival. This additional variable 
precludes the unambiguous determination of azimuthal angle 
from the single time difference equation. 

Furthermore, quadrantal ambiguities exist because 
of the periodic nature of the cosine function. For example, 
azimuth arrival angles of ^5°, 135°, 225°, and 315° would 
produce equal magnitude time differences between adjacent 
antenna responses. Fortunately, some differentiation 
between these responses is possible by knowing the algebraic 



sign of the time differences. This differentiation can be 
accomplished by a system using the phase coding technique 
because the chronological order in which the individual 
antennas respond to the incident wavefront will uniquely 
determine the composition of the phase coded sequences. 

For example, two incident wavefronts whose azimuthal arrival 
angles differ by l80° will create distinctively different 
phase coded sequences. Correlation filter pairs which are 
matched to either one of these sequences will not be matched 
to the other. 

This can be demonstrated using the example illustrated 
in Figures 2-8 and 2-11, with the rectangular pulse approach- 
ing from 180° instead of 0°. Assuming the same phase coding 
of the antenna signals, the combined antenna outputs for 
the 180° AOA signal will be as shown in Figure 2-lk . When 
those signals are input to the autocorrelation and cross- 
correlation filters (which remain matched for 0° AOA signals) 
the output waveforms will be as shovm in Figure 2-15. Com- 
bining the respective autocorrelation and crosscorrelation 
outputs from the channel A and channel B filters, the final 
output waveforms will be as shown in Figure 2-l6. These 
outputs indicate a severe mismatch between the coded signals 
and the correlation filters. Therefore, there is no sense 
ambiguity in the response of the system. In an actual RDF 
system utilizing this technique it would be possible to 
establish a bank of correlation filters which would respond 
(uniquely) to a full l80° range of azimuthal arrival angles. 
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However, there are still ambiguities created by 
azimuthal angles which produce equal (sign and magnitude) 
time differences. These ambiguities, and the ambiguities 
created by elevation angles of arrival must be resolved by 
an interferometer technique. A pair of linear arrays, 
situated in the ground plane and separated by a known fixed 
angle, will yield distinct, simultaneous time difference 
equations related to their respective baselines. A brief 
analysis of the two-array system, which is represented in 
Figure 2-17, will show that all spatial ambiguities can be 
resolved. 

The time difference equations which describe the 
(orthogonal) arrays in Figure 2-17 are: 

(At)^ = — cos cos a 
Id 1 

(At)„ = - cos cos a 

d C C. 

where the subscripted time shifts (At)^ and (At )2 are intro- 
duced by the corresponding linear arrays, and a represents 
the fixed elevation angle, which may take on any values 
between 0° and 90°. The azimuth angles 6^ and 02 are related 
by the orientation of the two linear arrays. For example, 
in the case of the orthogonal arrays in Figure 2-17; 



0 ^ = 



^7 



0 ^ + 270 °. 



These three equations will yield unambiguous solutions 
for the azimuth and elevation angles If the amplitude and 
sign of (At)^ and (At )2 can be determined. The uniqueness 
of these equations can be shown using the four Incident 
wavefronts represented In Figure 2-17 by the vectors I, II, 
III, and IV. An arbitrary elevation angle (a) Is assumed 
for all four wavefronts, and their azimuth angles are chosen 
as 30°, 150 °, 210°, and 330° relative to the axis of antenna 
number one. The time difference equations corresponding to 
these four wavefronts are as follows: 



Case I. (6^ = 30°, 

(At)^ = ^ cos 6 ^ 
(At)2 = cos 0 ^ 

Case II. (0^ = 150 °, 

(At)^ = ^ cos 0 ^ 

(At)., = — cos 0 ., 

2 c 2 

Case III: (0^ = 210°, 

(At), = — cos 0 , 

1 c 1 

(At)2 = - cos $2 



2 = 300°) 

cos a = +0.87 — cos a 
c 

cos a = + 0.50 — cos a 
02 = 60 °) 

cos a = -O.87 — cos a 
c 

cos a = + 0.50 — cos a 
c 

02 = 120 °) 

cos a = -O.87 — cos a 
c 

cos a = -O .50 ^ cos a 



i^8 



Case IV: (6^ = 330°, = 240°) 



(At )^ 


d 

= — cos 
c 


1 — 1 
CD 


cos a - + 0.87 


d 

— cos 
c 


(At )2 


d 

= — cos 
c 


^2 


cos a = -O. 5 O 


d 

— cos 
c 



Although the corresponding time shifts introduced 
by each of the v/avefronts have equal magnitudes, their 
algebraic signs are unique pair combinations. Therefore, 
differentiation between these (otherwise ambiguous) wave- 
fronts can be performed once the algebraic sign of each time 
shift is knov/n. That is, each pair of time difference 
values corresponds uniquely to a given azimuth angle and 
elevation angle combination. 

It has already been shown that the "sense" and 
therefore, the algebraic sign of the time shifts can be 
determined by the matched filter correlation process. 
Therefore, the two array system suggested here is theoret- 
ically capable of determining the parameters necessary to 
locate the incident wavefront in three dimensional space. 
However, the design of a specific mechanism to accomplish 
this is not attempted in this paper. The design, modeling, 
and evaluation of a physically realizable, three-dimensional 
RDF system utilizing phase coded linear arrays is considered 
a sufficient task for additional thesis work. 
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2 . Signal Pre-Selection 



As In any radio frequency receiving systems, the 
radio direction finder must be capable of electronically 
sorting through a range of operating frequencies In order 
to select particular signals of Interest. In most RDF 
systems currently In use, pre-selectlon and Identification 
of signals Is accomplished, either by a human operator or 
automatically, before the line of bearing Is obtained for 
that signal. The Idealized RDF system presented In this 
paper has thus far been limited to operation in a single- 
signal environment. However, in order to utilize the phase 
coding technique in a realizable system capable of operating 
In a dense RF environment, signal preselection and amplifi- 
cation functions must be incorporated into the signal 
processing portion of the system. 

The pre-selectlon method v;hich is used almost 
universally in communications receivers is the superhetero- 
dyne technique. That is, the desired signal frequency is 
converted by means of a non-linear mixing process to some 
Intermediate frequency (IF) at which narrowband high gain 
amplifiers may be used. In some cases the frequency con- 
verter circuit is preceded by one or more tuned RF amplifiers 
to provide coarse frequency selection, and to improve the 
receiver noise figure. However, the overall receiver 
selectivity is determined almost entirely by the IF circuits 
employed. 
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In order to apply the superheterodyne principle to 
the proposed RDF system, a decision must be made concerning 
where the frequency conversion should occur in the signal 
coding and decoding process. This decision will influence 
the type of circuits and components that would be required 
in the implementation of such a system, and in some cases 
may determine whether or not the system is operationally 
feasible. Basically, the conversion must occur at a stage 
of signal processing where it v/111 not affect the RDF opera- 
tion of the system. There are several locations in the 
signal path where frequency conversion and envelope detection 
might take place. 

The first possibility involves the phase coding 
and correlation of all input signals at the carrier fre- 
quency, with signal selection, conversion, and envelope 
detection at the output of the matched filters. The detected 
outputs of the matched filters would then be scanned for a 
unique crosscorrelation output. This approach has been 
utilized in the discussion and analysis of the phase coded 
system presented thus far in this paper. The sinusoidal 
analysis in the previous section has shown that phase coding 
and matched filter correlation of signals can be performed 
at the carrier frequency for all arrival angles. 

The primary disadvantage of this approach is that 
the signal processing is performed at very low signal levels. 
Accurate l80° phase inversion and amplitude balance of 
antenna outputs over a wide range of frequencies would 
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seem to be a very difficult requirement, considering that 
signal levels are typically in the microvolt range. These 
signal levels would be further reduced by the correlation 
filters, since most time delay devices, such as coaxial 
cable, electromechanical delay lines, and acoustic surface 
wave devices, are passive, and introduce signal losses. 

In addition, it is not clear what effect many 
different signals at various frequencies and arrival angles, 
will have on the phase coding and correlation processes. 

The possibility of signal interactions due to circuit 
non-linearities, which are always present in active devices, 
indicates that ambiguous results or spurious responses might 
be introduced by this v;lde-band approach. 

This approach does have the advantage hov/ever, that 
only one receiver is required. That is, the various antenna 
signals are eventually combined into a single channel which 
could then be input to a standard communications receiver. 

In this case the scanning process could be performed by 
electronically sampling the matched filter outputs and 
feeding the samples to the single receiver. 

A second possible configuration of the signal pro- 
cessor involves the frequency conversion of signals at each 
of the antenna outputs. With this procedure, the low level 
RP signals are converted to IP signals before the phase 
coding/decoding process occurs. In this case, the single- 
signal situation may again be assumed, with the coding/decoding 
being performed at the intermediate frequency Instead of 
the carrier frequency. 
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This approach would seem to offer several advantages 
over the previous wide-open front end approach. First of 
all, signal preselection before phase coding and correlation 
processing will eliminate the possibility of interaction 
between signals. Second, the noise figure of the system 
will be established by the individual array receivers, so 
the losses introduced by the correlation process vjIII have 
negligible effect upon the overall sensitivity of the 
system. Finally, since the phase coding and correlation 
process will be performed on high level, fixed frequency 
signals, more efficient narrowband circuits can be employed. 

The primary disadvantage of this system is that it 
requires the use of multiple receivers, which must be gain 
and phase matched over their total frequency range. However, 
it will be shown in the next chapter, that the use of large 
arrays is not necessary for good RDF resolution capability. 
Therefore, the use of short (2 or ^ element) array with 
this approach might result in a satisfactory trade-off 
between system performance and complexity. This approach, 
utilizing individual receivers at each antenna, and phase 
coding and correlation of signals at a fixed intermediate 
frequency, is employed in the system model evaluated (by 
computer simulation) in the following chapter. 
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Figure 2-1 , Phase Coded Input PiiLse and Matched Filter 
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Fisjure 2-2, Matched Filter Output Waveforms 
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Figure 2-3» Coinpleinentary Sequence Phase Coded Pulses and Matched 
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Figure 2-5, Complementary Sequence Phase Coded Pulses and Matched 
Crosscorrelation Filters 
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Figure 2-6, Crosscorrelation Filter Output V/aveforms 
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Figure 2-7, Four iilenent /mtenna Array and Induced Signals 
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Figure 2-8, Antenna Output Phase Coding Scheme 
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Figure 2-9, Antenna Output Waveforms 
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Figure 2-10, Combined Output Waveforms 
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A/C - Autocorrelation Output 
C/C - Crosscorrelation Output 



Figure 2-11, Four Element Correlation Filters 



64 






B (C/C) 



3 





t 




Figure 2-12, Correlation Filter Output Waveforms 
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Figure 2-1 3> Final Autocoirrelation and Crosscorrelation Outputs 
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Figure 2-1 4» Combined Antenna Outputs for 180° Arrival Angle 
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Figure 2-15, Correlation Filter Outputs for 180° Mismatch 
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Figure 2-1 6, Final Autocorrelation and Crosscorrelation Outputs 

For a 180*^ Mismatch 
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Figure 2-17, 



Two Array Interferometer 



III. SYSTEM MODELING 



In order to predict the effectiveness of the phase coded 
RDF technique in an actual signal environment, certain 
aspects of its performance should be evaluated. As in any 
RDF system, the most important parameters by which its 
usefulness can be judged are its angular resolution capabil- 
ity, and its performance in a noisy environment. It has 
already been shown that an idealized two array system is 
possible which can eliminate the gross bearing errors due 
to angular ambiguities. However, the basic performance 
parameters which will govern its precision are those Inherent 
to the individual linear arrays. Specifically, the bearing 
accuracy and resolution capability of the system is funda- 
mentally dependent on the precision of the signal processing 
technique used to determine time of arrival differences. 

Therefore a computer simulation of the signal encoding 
and decoding portion of the RDF system was employed to 
evaluate its performance under noisy signal conditions, and 
to determine the effects of varying array lengths and signal 
parameters. The program was originally developed by R. C. 
Todaro in his Master's thesis [1] to test the uniqueness of 
the correlation outputs for varying degrees of correlation 
filter mismatch. The original signal model used in the 
simulation was an Idallzed rectangular RF pulse of fixed 
duration, and a noise-free environment was assumed. The 
output from this simulation consisted of printer plots of 
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autocorrelation and crosscorrelation outputs for various 
arrival angles, and for varying degrees of correlation 
filter mismatch. 

The correlation results obtained from the original model 
Indicated that no Inherent ambiguities existed for arrival 
angles within any given quadrant of the array ground plane. 

In addition, the crosscorrelatlon plots obtained, using a 
l6-element array model. Indicated that the angular resolution 
of the system was limited only by the precision of the 
matched filter delay units. That Is, If delay lines could 
be physically realized so that their time delays corre- 
sponded to each Incremental degree of arrival angle, then 
the resolution of the system would be ±0.5°. 

A. MODEL DESCRIPTION 

For the purpose of this paper, the program was modified 
to produce continuous fine-grain autocorrelation and cross- 
correlatlon output plots from a CALCOMP plotter for varying 
noise levels, array lengths, and Input signal parameters. 

The system model Is Identical to the graphical Illustration 
used In the previous chapter. A single linear array was 
used, with the azimuthal angle of arrival measured from the 
baseline of the array. The signal model was modified as 
shown In Figure 3-1, to Investigate the effects of Input 
slgnal-to-nolse ratio on the correlation outputs. A 
constant noise level, corresponding to the average (or PIMS) 
level In an actual signal environment, was chosen. 
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PEAK SIGNAL LEVEL 



RMS NOISE LEVEL 



REFERENCE LEVEL 



Figure 3-1. Idealized Signal Model 



The signal-to-noise ratio was defined as the ratio of the 
peak signal voltage to the RMS noise voltage. 

The model parameters that are variable either by choice 
of specific program data cards, or by internal adjustments 
to program instructions, are as follov;s: 

1. Array length (2 , ^ , 8 , 10 , l6 , and l8 may be selected by 
data cards, larger arrays require re-dlmensloning of 
program storage). 

2. Arrival angle (selected by data cards) 

3. Input signal-to-noise ratio (selected by two program 
instructions). 

Choice of autocorrelation or crosscorrelation outputs 
(selected by data cards). 

5. Complementary sequences (any valid complementary code 
pair-may be loaded via data cards, the filter codes 
for autocorrelation or crosscorrelation processes are 
automatically derived by the program) . 
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6. Filter/Signal matching coefficient, D (the ratio of 
individual correlation filter delay to maximum 
propagation delay At, selected by data card). 

7. Input pulse v;idth (selected by program instruction). 

A detailed discussion of the program operation, format 

of input data, and an instruction listing are included in 
the appendix. 

B. MODEL OUTPUT RESULTS AND CONCLUSIONS 

The output plots obtained from this simulation are 
divided into three specific groups in order to show the 
effects of individual parameter changes on the overall 
operation of the model. These groups are outlined below, 
and are referenced to specific figure numbers , so that the 
related discussion and conclusions may be presented together 
at the end of this section. All pertinent parameters are 
listed in the heading of each of the plots. 

1. Group 1, Figures 3-2 through 3-^1 

These plots demonstrate the effect of varying the 
input signal-to-nolse ratio on autocorrelation and cross- 
correlation outputs. SNR takes on values of Infinity (no 
noise), 10:1, 5:1, and 2:1. 

Fixed parameters: Arrival Angle = 60.0° 

Array Length = l6 

Pulse Duration = 2.0 times maximum 

inter- antenna 
propagation time. 

For each combination of SNR and fixed parameters the matching 
coefficient (D) takes on values of 0.900, 0.950, 1.000, 1.050 
and 1.100. 
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2. Group 2 , Figures 3-^2 through 3-56 

These plots demonstrate the effect of varying array 
length on autocorrelation and crosscorrelation outputs. 

Array lengths of 10, 8 , and ^ elements are utilized. 

Fixed parameters: Arrival angle = 60.0° 

Pulse Duration = 2.0 times maximum 

inter-antenna 
propagation time. 

For each combination of array length and fixed parameters, 
the matching coefficient takes on values of 0 . 900 , 0 . 950 , 
1 . 00 , 1.050 and 1 . 100 . 

3 . Group 3 , Figures 3-57 through 3-80 

Effect of varying pulse width/antenna separation 
ratio on autocorrelation and crosscorrelation outputs. 

(Pulse durations of 20 and 200 times the inter-antenna 
propagation time are utilized) . 

Fixed parameters: Arrival Angle = 60.0° 

SNR = Infinite (no noise) 

4 . Discussion and Conclusions 

The plots comprising Group 1 are autocorrelation 
and crosscorrelation outputs for four decreasing values of 
signal-to-noise ratio. The first series, shown in Figures 
3-2 through 3-11 j is for an ideal noise-free signal environ- 
ment. The progression of the matching coefficient from 
0.900 to 1.100 in increments of 0.05 shows the expected 
result for the perfectly matched case (D = 1.000). The 
autocorrelation output is an exact replica of the input 
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pulse with a gain factor of 32, or tv:ice the array length. 

The crosscorrelation output is everyv;here zero for the 
matched case. 

As the slgnal-to-noise ratio begins to take on 
finite values, however, the crosscorrelation outputs no 
longer exhibit a zero value for the matched case. It is 
also evident that the autocorrelation outputs exhibit 
finite signal-to-noise ratios in the presence of input noise. 
A comparison of the matched autocorrelation outputs for 
input slgnal-to-nolse ratios of 10:1, 5:1, and 2:1 shov;s 
that the output signal-to-noise ratios are numerically 
equal to the input slgnal-to-noise ratios. 

The most important conclusion to be reached from 
this first group of plots, is that the crosscorrelation 
outputs still possess unique characteristics that distinguish 
them from the outputs of the mismatched filters, even in the 
presence of input noise. These characteristics are minimum 
amplitude and unipolar waveforms. The peak amplitudes of 
the crosscorrelation output waveforms show marked reductions 
as the perfectly matched case is reached. This characteris- 
tic is evident in the plots in Group 2 for array lengths of 
10, 8, and 4 as well. In general, for an input slgnal-to- 
nolse ratio of 5:1, all crosscorrelation outputs show an 
approximate 300^ reduction in peak output as the filter 
characteristic progresses from a mismatched condition to a 
perfectly matched condition (regardless of array length). 
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The second identifying feature of the crosscorrela- 
tion outputs is their polarity. For the mismatched case the 
outputs exhibit large amplitude noise peaks of both positive 
and negative polarity. V/hen the filters provide exact 
matching, however, the output is of a single polarity. For 
the l6-element array utilized in Group 1 these unipolar 
outputs are on the positive side. An interested quality of 
these unipolar outputs was noticed during the plotting of 
the various outputs for this paper. That is, the positive 
or negative polarity of these perfect match crosscorrelation 
outputs can be reversed by merely reversing the order of the 
complementary sequences used to phase code the input signals. 
While no attempt is made here to justify this result, it is 
mentioned because it might be of value in the design of 
recognizers for scanning the various correlation filter 
outputs . 

It should also be noted that the matched cross- 
correlation outputs do not possess unipolar qualities for 
all array lengths. For instance, the matched crosscorrela- 
tion output (with input noise) for an array length of 10, 
as shown in Figure 3-^^, is bipolar. It was assumed, there- 
fore, that this characteristic is dependent upon the chosen 
code length. Not enough data is available here to attempt 
any further explanation of this characteristic. 

The plots comprising Group 2 are 'autocorrelation and 
crosscorrelation outputs for array lengths of 10, 8, and 
4 elements. These plots were obtained in order to determine 
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the effect of array length on the resolution capability of 
the system. These plots, and those included in Group 1, can 
be used to derive an approximate maximum value of the 
resolution. That is, for matching coefficients of 0.950 
and 1 . 050 , which correspond to filters matched to arrival 
angles 1.7° above and below the true arrival angle, the 
correlation outputs are indicative of a mismatched condition. 
Therefore, the resolution of the system is at most 1.7°. 
Although finer resolutions are probable, and more precise 
limits could be assigned to this parameter by comparison of 
additional output plots, it is important to note that the 
resolution is not degraded by the shorter array lengths. 

The basic conclusion reached from these results was that a 
system utilizing short (2 or 4 element) arrays is technically 
feasible. It would seem that the only advantage to be de- 
rived from using longer array lengths is the increased signal 
processing gain in the autocorrelation mode. 

The plots comprising Group 3 are autocorrelation and 
crosscorrelation outputs for pulse widths normalized to 
20 and 200 times the maximum inter-antenna propagation time. 
This propagation time is equivalent to the time difference 
between antenna responses for a wavefront Incident at an 
arrival angle of zero degrees. It is obvious from these 
plots that the resolution of the system is degraded for 
signals of increased duration. The crosscorrelation outputs 
which were previously characteristic of the perfectly 
matched filter, now result for a range of filter mismatch. 
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For 20X pulse durations this output occurs over a range of 
filter mismatches from D = 1.000 (60° match) up to D = I. 050 , 
which corresponds to an arrival angle of 58.3°. For the 200X 
pulse this range of ambiguous outputs runs from 60 °, for a 
normally perfect match,' up to a matching coefficient 
corresponding to 25.8°. This last series of output plots 
indicates that the proposed system would not be effective 
against signals with long duration modulation events. 
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AUTOCORRELATION PLOT: AOA = 60,0°, FILTER DELAY = 0.900 TllffiS THE 

ANTENNA DELAY (FILTER IS MTCHED FOR 63.3°), ARRAY LENGTH = 16, 
SIGNALAOISE = INF., PULSE DURATION - 2.0 TIMES THE AJITEMA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE V/IDTH 100 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITSAlV, AMPLITUDE SCALE: 10 UNITSAlV. 



Figure 3-2 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0,950 TllffiS THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 61,7°), ARRAY LENGTH = 16, 
SIGNALA’OISE = INF,, PULSE DURATION =2.0 TllffiS THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORl^ALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORI'IALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITSAlV. AI-IPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-3 
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autocorrelation PLOT; AOA = 60.0°, FILTER DELAY = 1,000 TIMES THE 
ANTENNA DELAY {FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH = 16, 
SIGNAL/NOISE INF,, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 
















































































TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UIIETS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/1)IV. AMPLITUDE SCALE; 10 UNITS/DIV„ 



Figure 3-4 
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AUTOCORRELATION PLOT: AOA ^ 60,0°, FILTER DELAY 1,050 TII-IES THE 

ANTENNA DELAY (FILTER IS MTCHED FOR 58.3°), ARRAY LENGTH = 16, 
SIGNALAOISE = INF,, PULSE DURATION = 2,0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UlLETS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AlffLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. Al-IPLITUDE SCALE; 10 UNITS/t)IV. 



Figure 3-5 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY 1.100 TIMES THE 

ANTENNA DELAY {FILTER IS MATCIffiD FOR 56.6°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = INF., PULSE DURATION - 2.0 TllffiS TPIE ANTH1NA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE VDEDTH 100 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/DIV, Al-IPLITUDE SCALE: 10 UNITS/t»IV, 



Figure 3-6 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.900 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = INF., PULSE DURATION = 2.0 TIMES THE AIITENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE N0R14ALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UiriTS/DIVo AMPLITUDE SCALE: 5 UNITS/blV. 



Figure 3-7 
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CROSSCORRELATION PLOT; AOA = 60,0°, FILTER DELAY = 0,950 TIMES THE 
ANTEIWA DELAY (FILTER IS MATCHED FOR 6l.7°) ARRAY LENGTH = 16 , 
SIGNAL/NOISE - INF,, PULSE DURATION = 2,0 TD-IES THE ANTEIRIA DELAY, 




TIME (HORIZ) UNITS ARE NORM/HJZED TO IIH'UT PULSE WDTH 100 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TIME SCALE; 200 UNITS/DIV, AI-IPLITUDE SCALE; 5 UOTTS/blV, 



Figure 3-8 



CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.000 TUffiS THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH = 16, 
SIGNALAoISE --- INF., PULSE DURATION 2.0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE IflDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/t>IV. AMPLITUDE SCALE: 5 UNITSAlV. 

Figxire 3-9 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1 .050 TllffiS THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 58.3°), /iRRAY LENGTH r-. 16, 
SIGNALAoISE = IHF., PULSE DURATION = 2.0 TII4ES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITSAlV. AMPLITUDE SCALE; 5 UNITSAlV, 



Figure 3-10 



CROSSCORRELATION PLOT; AOA = 60.0°, FILTER DELAY = 1,100 TllfflS THE 
ANTENNA DELAY (FILTER IS MATCHED FOR 56,6°), ARRAY LENGTH = 16, 
SIGNAL/^OISE = INF,, PULSE DURATION = 2,0 TIl-IES THE AJJTENNA DELAY, 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE VJIDTH 100 UNITS 

EQUAL Olffi PULSE, AMPLITUDE UNITS ARE NORl-lALIZED TO PULSE HEIGHT, 
TIME SCALE; 200 UNITS/blV, AMPLITUDE SCALE; 5 UNITS/DIV, 



Figure 3-11 
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autocorrelation PLOT: AOA 60.0°, FILTER DELAY = 0.900 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 16, 
SIGNAL/^OISE = 10.0, PULSE DURATION =2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE V/IDTH 100 IHIITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/t)IV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-12 
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AUTOCORRELATION PLOT: AOA = 60 . 0 °, FILTER DELAY = 0.950 TIMES THE 

antenna DELAY (FILTER IS MATCHED FOR 61.7°), ARRAY LENGTH = 16, 
SIGNALAOISE = 10 . 0 , PULSE DURATION = 2,0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE IflDTH 100 UNITS 

EQUAL Olffi PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/blV. /i]»IPLITUDE SCALE: 10 UNITS/1)IV. 



Figure 3-13 
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autocorrelation PLOT: AOA = 60.0°, FILTER DELAY = 1.000 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH - 16, 
SIC2JAL/N0ISE - 10.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE VJIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 IRIITS/DIV. 



Figure 3-14 
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AUTOCORflELATION PLOT: AOA = 60.0°, FILTER DELAY = 1,050 TllffiS THE 

antenna delay (filter is matched for 58.3°), ARRAY LENGTH = 16, 
signal/noise = 10.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY, 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNlTS/DIV. AMPLITUDE SCALE: 10 UNITS/dIV. 



Figtire 3-15 
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autocorrelation PLOT; AOA = 60.0°, FILTER DELAY = 1.100 TDffiS THE 
antenna delay (FILTER IS MTCHED FOR $6.6°), ARRAY LENGTH = 16, 
signal/noise = 10.0, PULSE DURATION = 2.0 TIl-IES THE ANTENNA DEUY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE VfIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UiaTS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AMPLITUDE SCALE; 10 UNITS/tllV. 



Figiire 3-16 



CROSSCORRELATION PLOT: AOA 60.0°, FILTER DELAY = 0.900 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 10.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY, 




TIME (HORIZ) UIIITS ARE NORMALIZED TO INPUT PULSE Ii/IDTH 100 UNITS 

EQUAL ONE PULSE. Al^LITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/blV. AMPLITUDE SCALE: 5 UNITS/DIV. 



Figure 3-17 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY 0.950 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 6l.7°), ARRAY LENGTH = 16, 
SIOIALAOISE = 10.0, PULSE DURATION = 2.0 TB4ES THE ANTENNA DELAY. 




HME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE V/IDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE TOUTS ARE NORMALIZED TO PULSE HEIGHT. 
nME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 5 UNITS/DIV. 



Figure 3-18 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY - 1.000 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 60,0°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 10.0, PULSE DURATION = 2.0 TIMES THE Al^TENNA DELAY. 




TIME (HORIZ) UNITS ARE NORI^LIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AMPLITUDE SCALE; 5 UNITS/DIV. 



Figure 3-19 



CROSSCORRELATION PLOT: AOA = 60,0°, FILTER DELAY = 1.050 TIMES THE 

antenna delay (filter is MATCHFD for 58.3°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 10.0, PULSE DURATION = 2.0 TIMES THE ANIEITNA DELAY. 




HME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE LTEDTH 100 UNITS 

equal ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/lDIV. AMPLITUDE SCALE: 5 UNITS/tllV. 



Figiure 3-20 



CROSSCORRELATION PLOT; AOA = 60.0°, FILTER DELAY = 1.100 TIJffiS THE 
ANTENNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 10.0, PULSE DURATION =2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AMPLITUDE SCALE: 5 UOTTS/tolV. 



Figure 3-21 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.900 TIMES THE 

antenna DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 16 , 
signal/noise = 5.0, PULSE DURATION = 2.0 TIMES THE AI^TEMA DELAY, 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE ;flDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/t)IV, 



Figure 3-22 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY - 0.950 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 61,7°) ARMY LENGTH = 16, 
SIGNALAOISE = 5.0, PULSE DUR/vTION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORl-IALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. Al-IPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/1)IV. Al-IPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-23 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY -- 1.000 TIMES THE 

ANTEIWA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH = 16, 
SICaiAL/^OISE = 5.0, PULSE DURATION = 2.0 TIMES TJIE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE N0RI4ALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/blV. AMPLITUDE SCALE: 10 UNITS/blV. 



Figiire 3-24 



AUTOCORRELATION PLOT: AOA - 60,0°, FILTER DELAY -- 1.050 TIMES THE 

antenna DELAY (FILTER IS MiATCHED FOR 58.3°), ARRAY LEtIGTH = 16, 
SIGNALAoISE = 5.0, PULSE DURATION = 2.0 TIMES THE ANTEMA DELAY. 




TIME (HORIZ) UILETS ARE NORMALIZED TO INPUT PULSE WDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/tolV. 



Figure 3-25 
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AUTOCORRELATION PLOT: AOA = 60,0°, FILTER DELAY - 1.100 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 16, 
SIGNAL/NOISE r= 5 . 0 , PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/t)IV. 



Figure 3-26 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.900 TIMES THE 

antenna delay (filter is matched for 63 . 3 °), ARRAY LENGTH = 16, 
SIGNAL/4J0ISE = 5.0, PULSE DURATION = 2.0 TIl-IES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE IflDTH 100 UNITS 

EQUAL ONE PULSE. Al-IPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/blV. AMPLITUDE SCALE; 5 UNITS/blV. 



Figure 3-27 
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CROSSCORRELATION PLOT; AOA = 60.0°, FILTER DELAY = 0.950 TB'IES THE 
AI^TENNA delay (FILTER IS MATCHED FOR 6l.7°), ARRAY LENGTH = 16, 
SIGNALAOISE = 5.0, PULSE DURATION = 2.0 TIMES THE AfITENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO IIIPUT PULSE WIDTH 5 UNITS 

EQUAL ONE PULSE. AI-IPLITUDE UNITS ARE NORl'IALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. Al^LITUDE SCALE; 5 UNITSAlV. 



Figure 3-28 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.000 THIES THE 

ANTEOTIA delay (filter is matched for 60.0°), ARRAY LENGTH = 16, 
SIGNAL/^IOISE = 5.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITSA>IV. AMPLITUDE SCALE: 5 UNITS/DIV, 



Figure 3-29 
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CROSSCORRELATION PLOT; AOA = 60,0°, FILTER DELAY = 1,050 TIMES THE 
ANTENNA DELAY (FILTER IS MATCHED FOR 58.3°), ARRAY LENGTH = 16 , 
SIGNAL/^OISE = 5.0, PULSE DURATION = 2,0 TIMES THE Al'fTENIW DELAY. 
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TIME (HORIZ) UNITS ARE NORI^ALIZED TO INPUT PULSE V/IDTH 100 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TDffi SCALE: 200 UNITS/DIV, A1-IPLITUDE SCALE: 5 UNITS/DIV. 



Figure 3-30 



CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.100 TllffiS THE 

ANTENNA DELAY (FILTER IS MATCFIED FOR 56.6°), /iRRAY LENGTH = 16, 
SIGNAL/NOISE = 5.0, PULSE DURATION = 2.0 TI14ES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORl-lALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 5 UNITS/t»IV. 



Fig\u?e 3-31 
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AUTOCORRELATION PLOT: AOA = 60.0°, nLTER DELAY 0.900 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH 16, 
SIGNALAOISE = 2,0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




HME (HORIZ) UNITS ARE NORMALIZED TO IITPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-32 



AUTOCORRELATION PLOT: AOA = 60.0°, FIL'IER DELAY = 0.950 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 61.7°), ARR/.Y LENGTH == 16, 
SIGNAL/NOISE - 2.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORI-IALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/DIV, AMPLITUDE SCALE: 10 UNITS/DIV, 



Figure 3-33 
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AUTOCORRELATION PLOT; AOA = 60.0°, FILTER DELAY = 1.000 TIMES THE 
ANTENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH 16, 
SIGNAL/NOISE =2.0, PULSE DURATION =2.0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE VJIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORl-lALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/i)IV. AMPLITUDE SCALE; 10 UNITS/DIV. 



Figure 3-34 
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AUTOCORRELATION PLOT; AOA ^ 60.0°, FILTER DELAY = 1.050 TIMES THE 
ANTENNA DELAY (FILTER IS MATCHED FOR 58.3°), ARRAY LENGTH = 16, 
signal/noise = 2.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/1)IV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-35 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.100 TIMES TJIE 

ANTENNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 2.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE TODTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV, AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-36 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0,900 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 16, 
SIGNALAOISE = 2.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORI-IALIZED TO INPUT PULSE \flDTH 100 UfJITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORI-IALIZED TO PULSE HEIGHT. 
TIl-lE SCALE: 200 UNITS/1)IV. AMPLITUDE SCALE: 5 UNITS/l)IVo 



Figure 3-37 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.950 TllffiS THE 

ANTEIWIA DELAY (FILTER IS MATCHED FOR 61.7°), ARRAY LH1GTH = 16, 
SIGNAL/^IOISE = 2.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORflALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMAI.IZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/blV. AMPLITUDE SCALE: 5 UNITS/blV. 



Figure 3-38 
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CROSSCORRELATION PLOT: AOA = 60*0°, FILTER DELAY = 1,000 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 2.0, PULSE DURATION = 2.0 TIMES THE AI'ITENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 5 UNITS/t)IV. 



Figvire 3-39 



CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.050 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 58.3°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 2.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE ’NIDTH 100 UOTTS 

EQUAL ONE PULSE. AMPLITUDE UIJITS ARE NORITALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/tUV. AI-IPLITUDE SCALE: 5 UNITS/t)IV. 



Figure 3-40 
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CROSSCORRELATION PLOT: AOA 60.0°, FILTER DELAY = 1.100 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = 2.0, PULSE DURATION = 2.0 TIMES THE ANTEMA DELAY, 




TIME (HORIZ) UNITS ARE N0R1-IALIZED TO INPUT PULSE V/IDTH 100 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 Ul^ITS/DIV, AMPLITUDE SCALE; 5 UNITS/blV, 



Figure 3-41 
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CROSSCORRELAnON PLOT; AOA = 60.0°, FILTER DELAY = 0,900 TI1>ES TliE 
ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 10, 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV, AMPLITUDE SCALE: 5 UiaTS/DIV, 



Figure 3-42 



CROSSCORRELATION PLOT: AOA = 60.0°, ITiLTER DELAY = 0.950 TDffiS THE 

ANTEMA DELAY (FILTER IS MATCHED FOR 6l.7°), ARRAY LENGTH = 10, 
SIGNAL/NOISE ^ 5 . 0 , PULSE DURATION = 2.0 TIMES THE ANTE^mA DELAY. 




EQUAL ONE PULSE. AMPLITUDE IRIITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/blV. AMPLITUDE SCALE: 5 UNITS/1)IV. 



Figure 3-43 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.000 TIMES THE 

AlWENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH 10, 
SIGNAL/HOISE = 5.0, PULSE DURATION = 2.0 TIMES THE ANTETTNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO Il^PUT PULSE VttDTH 200 UNITS 

EQUAL ONE PULSE. Al^lPLITUDE IHTITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. Al-IPLITUDE SCALE: 5 UNITS/l)IVo 



Figure 3-44 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1 .050 TII-IES THE 

ANTEIRIA DELAY (FILTER IS MATCHED FOR 58.3°), ARRAY LENGTH = 10, 
SIGNAL/NOISE = 5.0, PULSE DURATION = 2.0 TIMiES THE AlITENNA DELAY. 




TIME (HORIZ) UNITS ARE NORt'IALIZED TO INPUT PULSE V/IDTH 200 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORI-IALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 5 UNITS/tHV, 



Figure 3-45 



CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.100 TIMES THE 

ANTENNA DELAY (HLTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 10, 
SIGI^ALAOISE = 5.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL Olffi PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/t)IV. AMPLITUDE SCALE: 5 UNITS/DIV. 



Figure 3-46 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.900 TI14ES THE 

ANTENNA DELAY (FILTER IS MTCHED FOR 63.3°), ARRAY LENGTH = 8, 
SIGNAL/NOISE = 5.0, PULSE DURATION = 2.0 TIMES THE AI^TENNA DELAY. 




TIME (HORIZ) UNITS ARE NORI^IZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE<> AMPLITUDE UNITS ARE NORILiLIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/DIV, Al-IPLITUDE SCALE: 2 UNITS/DIV. 



Figure 3-47 
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CROSSCORRELATION PLOT: AOA - 60,0°, FILTER DELAY = 0.950 TIMES THE 

ANTENNA DELAY (FILTER IS MTCHED FOR 61,7°), ARRAY LENGTH = 8, 
SIGNAL/NOISE ^ 5,0, PULSE DURATION = 2,0 TIl-ES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE lOTTH 200 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORNIALIZED TO PULSE HEIGHT p 
TIME SCALE: 200 UNITS/t)IV, AMPLITUDE SCALE: 2 UNITS/DIV, 



Figure 3-48 



CROSSCORRELATION PLOT; AOA = 60.0°, FILTER DELAY = 1.000 TIMES THE 
ANTENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH = 8, 
SIGNAL/NOISE = 5.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE. Al-IPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/t)IV. AMPLITUDE SCALE: 2 UNITS/DIV. 



Figure 3-49 
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CROSSCORRELATION PLOT: AOA = 60 . 0 °, FILTER DELAY = 1,050 THffiS THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 58.3°), ARRAY LENGTH = 8, 
SIGNAL/NOISE = 5.0, PULSE DURATION = 2,0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORIiALIZED TO PULSE HEIGHTo 
TIME SCALE: 200 UNITS/DIV, AMPLITUDE SCALE: 2 UNITS/blV, 



Figure 3-50 
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CROSSCORRELATION PLOT: AOA 60.0°, FILTER DELAY = 1,100 TI^^ES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 56,6°) , ARRAY LENGTH = 8, 
SIGNAL/NOISE =5.0, PULSE DURATION =2.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UHITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 2 UNITS/DIV. 



Figure 3-51 
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CROSSCORRELATION PLOT: AOA = 60,0°, FILTER DELAY ^ 0.900 TllffiS THE 

ANTENNA DELAY (FILTER IS MTCHED FOR 63.3°), ARRAY LENGTH = 4, 
SIGJJAL/NOISE = 5.0, PULSE DURATION = 2.0 TIl-IES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO IIJPUT PULSE raDTH 200 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORl'IALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 1 UNIT/DIV. 



Figure 3-52 
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CROSSCORRELATION PLOT: AOA 60.0°, FILTER DELAY - 0.950 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 61,7°), ARRAY LENGTH = 4, 
SIGNALAOISE = 5.0, PULSE DURATION = 2.0 TIMES THE AJfTENNA DELAY, 
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time (HORIZ) units are NORMALIZED TO INPUT PULSE V/IDTH 200 UI^TS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORI-IALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/DIV, AMPLITUDE SCALE: 1 UNIT/t)IV, 



Figure 3-53 
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CROSSCORRELATION PLOT; AOA = 60.0°, FILTER DELAY =•- 1,000 TB4ES THE 
ANTEI'INA DELAY (FILTER IS MATCHED FOR 60,0°), /JIRAY LENGTH = 4, 
SIGI'JAL/NOISE = 5.0, PULSE DURATION = 2,0 TllffiS THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE, AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TIl'lE SCALE: 200 UNITS/t)IV. AMPLITUDE SCALE; 1 UNIT/DIV. 



Figure 3-54 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.050 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 58.3°), ARRAY LENGTH = 4 , 
SIGNAL^OISE 5.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE. A1-5PLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIl-IE SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 1 UNIT/DIV. 



Figure 3-55 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.100 TIMES THE 

ANTEITNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 4 » 
SIGNAL/i'IOISE = 5.0, PULSE DURATION = 2.0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORtlALIZED TO INPUT PULSE WIDTH 200 UNITS 

EQUAL ONE PULSE. AI^LITUDE UNITS ARE NORI-IALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/1)IV. AI'IPLITUDE SCALE: 1 UNIT/1)IV. 



Figure 3-56 
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AUTOCORRELATION PLOT: AOA = 60 . 0 °, FILTER DELAY = O. 6 OO TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 72.6°), ARRAY LEITOTH = I 6 , 
signal/noise = INF., PULSE DURATION =20.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 4 OO UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-57 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.800 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 66.4°), ARRAY LEl'IGTH = 16, 
SIGNALAOISE = INF., PULSE DURATION = 20.0 TIMES THE ANTENNA DELAY, 
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TIME (HORIZ) UNITS ARE NORf-lALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL ONE PULSE. AI4PLITUDE UlTITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UlLETS/DIV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-58 



AUTOCORRELATION PLOT: AOA = 60,0°, FILTER DELAY = 0,900 TII4ES THE 

ANTENl^A DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH r-. 16, 
SIGNAL/lJOISE = INF., PULSE DURATION = 20.0 TIMES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE N0R2-1ALIZED TO INPUT PULSE VJIDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORl-IALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. Al'SPLITUDE SCALE: 10 UNITS/DIV„ 



Figure 3-59 
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autocorrelation PLOT: AOA = 60 . 0 °, FILTER DELAY = 0.950 TIMS TM 

ANTENNA DELAY (FILTER IS MTCHED FOR 61.7°), ARRAY LEI'TGTH = 16, 
signal/noise = INF., PULSE DURATION = 20.0 TD4ES THE AffTENNA DELAY 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UlOTS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AI-IPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-60 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.000 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 60.0°), /lRRAY LENGTH - 16, 
SIC^AL/NOISE = INF., PULSE DURATION = 20.0 TllffiS THE ANTENNA DELAY, 
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time (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE lODTH 4 OO UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHTc 
HME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-61 
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AUTOCORRELATION PLOT: AOA = 60,0°, FILTER DELAY = 1.050 TIMES THE 

ANTENNA DELAY (PILTER IS MATCHFD FOR 58.3°), AIUIAY LENGTH = 16, 
SIGNAL/NOISE = INF., PULSE DURATION =20.0 TIMES THE ANTENNA DELAY, 
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TDffi (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS i\RE NORMALIZED TO PULSE HEIGHT. 
HME SCALE: 200 UNITS/DIV. AI^LITUDE SCALE: 10 UNITS/DIV. 



FigTxre 3-62 
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AUTOCORRI'IATION PLOT: AOA = 60,0°, FILTER DELAY = 1,100 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = INF., PULSE DURATION = 20.0 TIMES THE ANTHWA DELAY, 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE ICLDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/t)IV, AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-63 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.200 HMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 53.1°), ARRAY LENGTH =•: 16, 
SIGNAL/NOISE = INF., PULSE DURATION = 20.0 TIMES THE ANTEiniA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE IflDTH — 400 UNITS 
EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-64 
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AUTOCORRELATION PLOT: AOA = 60.0^^, FILTER DELAY = 1 . 4 OO TBIES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 45.6°), ARRAY LENGTH = I 6 , 
SIGNAL/NOISE = INF., PULSE DURATION = 20.0 TBIES THE ANTENNA DELAY. 




TIME (HORIZ) UNITS ARE NOR>LALIZED TO INPUT PULSE V/IDTH 40 O UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE N0R1>IALIZED TO PULSE HEIGHT. 
TUffi SCALE: 200 UNITS/blV, Al-IPLITUDE SCALE: 10 UNITS/DIV. 



Figure 3-65 



CROSSCORRELATION PLOT: AOA 60.0°, FILTER DELAY = 0.900 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 63.3°), ARRAY LENGTH = 16, 
SIGNAL/iWISE = INF., PULSE DURATION = 20.0 TIMES THE AITIENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL ONE PULSE, AI-IPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AI'IPLITUDE SCALE: 5 UNITS/t)IV. 



Figure 3-66 
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CROSSCORRELATION PLOT: AOA = 60,0°, FILTER DELAY = 0.950 TIKES THE 

AI^TEIRIA DELAY (FILTER IS MATCHED FOR 61.7®), ARRAY LENGTH 16, 
SIGNAL/HOISE - I NT,, PULSE DURATION = 20.0 TIMES THE ANTENNA DELAY. 




time (HORIZ) units are NORI'IALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORf-IALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AI-IPLITUDE SCALE; 5 UlMTS/llIV. 



Figure 3-67 
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CROSSCORRELATION PLOT: AOA 60.0°, FILTER DEUI = 1.000 TIMES THE 

AOTENNA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH 
SIGNAL/NOISE = INF., PULSE DURATION 20.0 TIMES THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE raOTH 400 UNITS 

EQUAL ONE PULSE. Al’lPLITUDE UNITS ARE NORMALIZED TO PUI.SE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. /J4PLITUDE SCALE: 5 U■NITS/DI\^ 



Figiire 3-68 
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CROSSCORRELATION PLOT: AOA = 60 , 0 °, FILTER DELAY = 1,050 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 58,3°), ARRAY LENGTH =16, 
signal/noise ■- INF,, PULSE DURATION = 20,0 TIMES THE ANTENNA DELAY, 
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TIME (HORIZ) UNITS ARE N0R1>IALIZED TO INPUT PULSE V/IDTH 40 O UI'flTS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 5 UNITS/DIVo 



Figure 3-69 



CROSSCORRELATION PLOT; AOA - 60,0°, FILTER DELAY = 1,100 TIMES TPIE 
AlPTErPNA DELAY (FILTER IS MATCHED FOR 56.6°), ARRAY LENGTH r. 16, 
SIGNAL/NOISE =•- INF., PULSE DURATION = 20.0 TllffiS TPIE MTENNA DELAY. 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE V/IDTH 400 UNITS 

EQUAL ONE PULSE, Al-IPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/1)IV. AMPLITUDE SCALE; 5 UNITS/DIV. 



Figure 3-70 
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autocorrelation PLOT: AOA = 60.0°, FILTER DELAY = 0.200 TIl'ES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 84.3°), ARRAY LENGTH = 16, 
SIGNAL/lTOISE = INF,, PULSE DURATION = 200.0 TllffiS THE Al'ITENNA DELAY 
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TIME (HORIZ) UI'IITS ARE NORJ-IALIZED TO INPUT PULSE VJIDTH 400 UNITS 

EQUAL ONE PULSE. AMPLIITIDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/blV. 



Figure 3-71 
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AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.600 TllffiS THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 72.6°), ARRAY LENGTH = 16, 
SIGNALA'OISE ^ INF., PULSE DURATION = 200.0 TIMES THE ANTEMA DELAY, 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/t)IV. AMPLITUDE SCALE: 10 UNITS/t)IV. 



Figure 3-72 
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AUTOCORRELATION PLOT; AOA - 60.0°, FILTER DELAY = 1.000 TIMES THE 
ANTENNA DELAY (FILTER IS MATCHED FOR OO.o'^), ARRAY LENGTH = 16, 
SICa^AL/NOISE - INF., PULSE DURATION = 200.0 Tll-ffiS THE AI'ITENNA DELAY, 




TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL 01® PULSE, AMPLITUDE UNITS ARE NORI-IALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/t)IV. AITLITUDE SCALE; 10 UNITS/blV, 



Figure 3-73 



AUTOCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.400 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 45.6°), ARRAY LENGTH - 16 , 
SIGNAL/NOISE = INF,, PULSE DURATION - 200.0 TII-IES THE ANTEMJA DELAY 
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TIME (HORIZ) UNITS ARE NORI-IALIZED TO HH^UT PULSE WIDTH A-00 UNITS 

EQUAL ONE PULSE. AMPLITUDE DiLLTS ARE N0R14ALIZED TO PULSE HEIGHT, 
TIME SCALE: 200 UNITS/DlY. AMPLITUDE SCALE: 10 UNITS/t)IV, 



Figiire 3-74 
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AUTOCORRELATION PLOT: AOA = 60,0°, FILTER DELAY = 1.800 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 25.8°), ARRAY LENGTH = 16, 
SIGNAL/iroiSE INF,, PULSE DURATION = 200.0 TIMES THE ANTEimA DELAY 
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TIME (HORIZ) UNITS ARE NORJ-IALIZED TO INPUT PULSE WIDTH 400 UNITS 

EQUAL ONE PULSE, Al'IPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UillTS/DIV. AMPLITUDE SCALE; 10 UNITS/DIV, 



Figure 3-75 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 0.200 TllffiS THE 

ANTEMA DELAY (FILTER IS FiATCHED FOR 84.3°), ARRAY LENGTH 16, 
SlffilAL/NOISE = INF., PULSE DURATION = 200.0 TIMES THE ANTEtfllA DELAY. 
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THUS (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE IHDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UlIITS/blV. Al^IPLITUDE SCALE; 10 UNITS/DIV. 



Figure 3-76 



CROSSCORRELATIOK PLOT; AOA = 60.0°, FILTER DEL.'lY - 0.600 TIMES TliE 
ANTENNA DELAY (FILTER IS MATCHED FOR 72.6°), ARIUY LENGTH =16, 
SIGNAL/lJOISE = INF,, PULSE DURATION = 200.0 TII«S THE ANTENNA DELAY. 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE \flDTH 400 UNITS 

EQUAL ONE PULSE, A1>1PLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/tHV, A14PLITUDE SCALE; 10 UNITS/t)IV. 



Figure 3-77 
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CROSSCORRELATION PLOT: AOA ^ 60.0°, FILTER DEL/iY = 1.000 TllffiS THE 

ANTE^WA DELAY (FILTER IS MATCHED FOR 60.0°), ARRAY LENGTH = 16, 
SIGNAL/NOISE = INF., PULSE DURATION = 200.0 TIMES THE ANTENNA DELAY 
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TIME (HORIZ) UNITS ARE NORMALIZED TO INPUT PULSE VEEDTH 400 UNITS 

EQUAL Olffi PULSE. Al'lPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGIiT. 
TIME SCALE; 200 UNITS/DIV. AMPLITUDE SCALE: 10 UNITS/t)IV. 

Figure 3-78 
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CROSSCORRELATION PLOT: AOA = 60.0°, FILTER DELAY = 1.400 TIMES THE 

ANTENNA DELAY (FILTER IS MATCHED FOR 45.6°), ARRAY LENGTH = 16, 
SIGNALAOISE =: INF., PULSE DURATION = 200.0 TIJ4ES THE ANTENNA DELAY, 
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TIME (HORIZ) UNITS ARE NOPIIALIZED TO INPUT PULSE klDTH 400 UNITS 

EQUAL ONE PULSE. AMPLITUDE UNITS ARE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE: 200 UNITS/DIV. AMPLITUDE SCALE: 10 UIMTS/lDIV. 



Figure 3-79 
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CROSSCORRELATION PLOT; AOA = 60.0°, FILTER DELAY 1.800 TIMES THE 
ANTENNA DELAY (FILTER IS MTCIiED FOR 25.8°), ARRAY LENGTH = 16, 
SIGIIALAOISE = INF., PULSE DURATION = 200.0 TIMES THE ANTEJRU DELAY, 
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TIME (HORIZ) UNITS ARE NORl'lALIZED TO INPUT PULSE WIDTH 400 UOTTS 

EQUAL ONE PULSE. AJ'IPLITUDE UNITS ilRE NORMALIZED TO PULSE HEIGHT. 
TIME SCALE; 200 UNITS/DIV. AMPLITUDE SCALE: 10 Ul-IITS/DIV. 



Figure 3-80 
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IV, SUIMARY AND 'RECOMMENDATIONS 



A. SUMMARY 

The phase coding radio direction finder using comple- 
mentary sequence pairs has been evaluated with the objective 
of applying the basic technique to a realizable model of the 
system. A sinusoidal analysis of the carrier phase relation- 
ships which exist In the various portions of the signal 
process has shown that the unique autocorrelation and cross- 
correlation outputs will be produced for any arrival angle. 

If the system Is matched for that angle. 

It was shov'/n that the Inherent ambiguities of a single 
array system., caused by elevation angle and adjacent quadrant 
azimuth components, can be resolved by using a common crossed- 
array Interferometer technique. This technique requires that 
a set of correlation filters corresponding to a l80° range 
of arrival angles be available for each array. 

A computer model of the system was used to evaluate the 
performance of the signal processor under varying Input noise 
conditions, and for varying array lengths and signal duration 
parameters. The conclusions reached from this analysis 
Indicate that the angular resolution of the system will not 
be seriously affected by the array length or Input slgnal-to- 
nolse ratio. That Is, for most Input signals, the resolving 
capability of the system Is dependent only upon the number 
and Incremental differences between the time delay segments 
of the correlation filters used. However, It was also 



159 



determined that the duration of the signal model used for 
this analysis had a very significant effect upon the resolu- 
tion. For signal pulse v;ldths that were large compared with 
the time of arrival difference between antennas, the resolu- 
tion of the system V'/as seriously degraded. As a result, it 
was concluded that the system would not be effective against 
signals with this characteristic. 

B. RECOMMENDATIONS FOR FURTHER INVESTIGATION 

Since the analysis of the proposed system v;as based 
entirely upon the relative time of arrival of distinct 
modulation events, and since this was shovm to place a 
serious limitation on the type of modulation signals that 
can be processed by the system, an alternate technique 
should be pursued. That is, the utilization of carrier 
phase differences, resulting from the propagation delay 
between antennas, might be a more acceptable parameter to 
use rather than time difference. 

A second recoirimendat ion v:ould be the construction of a 
simple scale model, utilizing one or two small arrays in 
order to evaluate the operation of the system in a multi- 
signal environment. 
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APPENDIX 



COMPUTER PROGRAM 

The computer simulation is a digital program written in 
FORTRAN IV level G. It was processed by an IBM 36O/67 
system and the built in functions are those normally 
available on this system. 

The program requires a data deck input for the specific 
complementary code used in addition to the parameters 
required to process the code. The following words are 
peculiar to this program. 



CODEl: 


First sequence in a complementary code 


C 0 DE 2 : 


Second sequence in a complementary code 


SIGN: 


Summation sign, used for addition of correlator 
outputs 


DELDEG: 


Incremental angle shift in degrees, used to 
obtain a new value for the angle 


D: 


Filter delay starting value 


START: 


The starting value in degrees of the incoming 
signal 


STOP : 


The final value in degrees of the incoming signal 


START-T: 


The starting time for the computation 


STOP-T: 


The completion time for the computation 


DSTEP: 


Incremental step for the parameter D, used to 
obtain a new value for the filter (correlator) 
delay 


NRANTS : 


The number of antennas in use 


SPACE: 


Spacing between antennas; antennas have uniform 
spacing 


DEBT: 


Computation time Increment 
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TIMES: Multiplying factor used in the subprograms In 

order to utilize Integer arithmetic and logic 

JJ: Control for suppressing data printout; do not 

print if JJ--2 

JJJ: Control for correlation function; if JJJ=1 then 

cross-correlate j if not then autocorrelate 

JJJJ: Control to suppress the final plot v/hich is a 

summary of maximum amplitude versus angle; if 
JJJJ=2 then do not plot 

The data deck must be assembled in a specific sequence 
in order to have the program operate properly. Additionally 
if the data is not known a card must still be loaded into 
the data deck in order not to invalidate the read instruc- 
tions. The data must be provided whether valid or not. 

The data cards are assembled as in Figure A-1. The first 
card in the data deck is CODEl and the last is JJJJ. All 
data cards are punched starting in column one. CODEl and 
C0DE2 are formatted floating point . 1 and each card can 
have up to eighteen entries. Each code is read from one 
card, if longer codes (greater than l8) are required then 
the read statement must be changed to read more than one 
card for each code. NRANTS, JJ, JJJ and JJJJ are formatted 
in integer 111 and all other data cards are formatted in 
floating point F11.6. 

The program can provide for a pulse signal input with 
the pulse duration equal to any multiple of the delay due 
to antenna separation. The program will automatically scale 
the built in pulse duration (W) and the antenna separation 
(SPACE) to accommodate any code length. Additionally the 



162 



computation time is automatically scaled in order to 
include all the correlation outputs greater than zero and 
also provide for the maximum number of computations per 
angle and filter delay setting. The maximum number of 
computation intervals for a fixed angle and delay is 2000. 

The codes are read into computer memory and all necessary 
operations are accomplished in order to provide the proper 
autocorrelation and crosscorrelation functions as a function 
of time, angle and filter delay. If data is provided which 
is beyond the array storage set aside for the computation 
then the program will attempt adjustments in the size of the 
data. After twenty attempts vjlth no success the program will 
terminate and print out: DATA EXCEEDS ARRAY STORAGE... 

JOB TERMINATED. 

The program v/ill input all parameters, compute other 
required parameters in accordance with the inputted param- 
eters, compute an output for a fixed time, angle and delay 
and then step the time variable to the next increment. After 
the time variable has been stepped across its range the 
filter matching coefficient (D) is Incremented to a new 
value, and the output is computed for all time values as the 
time variable is again swept through its range. After the 
matching coefficient has been swept through values from its 
initial value to a final value equally spaced above D=1.000, 
the angle variable may then be incremented to a new value. 

The outputs are then computed across the time range and the 
angle range. In this manner the outputs of the filters 
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(correlator outputs) are computed for all times of interest, 
all angles of interest and all desired variations of the 
filter delay lines. 

The outputs available from this program are printer 
plots and fine grain continuous plots of correlation output 
versus time for each and every value of angle and filter 
delay. The printer plots contain all pertinent parameters 
and scaling Information. 

A listing of important parameters and statements con- 
cerning the Interpretation of the various control commands 
are printed at the beginning of each computer run tp aid in 
analysis of the data and the plots. The maximum core size 
is I5OK and, due to the fine grain computation required by 
the program, execution time for most investigations is in 
the 15 to 20 minute region. 
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TYPICAL DATA DECK 



1.0 1.0 1 . 0 - 1. 0 1.0 1 . 0 - 1.0 1.0 1.0 1.0 1 . 0 - 1 . 0 - 1 . 0 - 1.0 1 . 0 - 1.0 

1.0 1.0 1 . 0 - 1.0 1.0 1 . 0 - 1.0 1 . 0 - 1 . 0 - 1 . 0 - 1.0 1.0 1.0 1 . 0 - 1.0 1.0 

1.000000 

3.000000 
0.990000 

60.000000 

60.000000 
0.000000 
0.000000 
16 
2 
2 
2 



CARD 1. . .FIRST SEQUENCE OF THE CODE PAIR (CODEl) 

CARD 2... SECOND SEQUENCE OF THE CODE PAIR (CODE2) 

CARD 3. . .SUMMATION SIGN (SIGN) 

CARD . .ANGLE STEP (DELDEG) 

CARD 5... filter DELAY STARTING VALUE (D) 

CARD '6. . .ANGLE START (START) 

CARD 7... angle stop (STOP) 

CARD 8. . .COMPUTATION TIME START (START-T) 

CARD 9... filter delay STEP (DSTEP) 

CARD 10.. NUMBER OF ANTENNAS (NRANTS) 

CARD 11.. CONTROL FOR DATA PRINT (JJ) 

CARD 12.. CONTROL FOR TYPE OF CORRELATION (JJJ) 

CARD 13.. CONTROL FOR OUTPUT VS ANGLE PLOT (JJJJ) 
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TC ENSURE THAT CORRELATION PROCESS CAN BE CCNPLETED IN 
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LCAC CODE 1 INTO THE DUMMY ARRAY 
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